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Abstract 

Some  efficient  new  classes  of  algorithms  for  compensating  for  fading  in  wireless  systems 
are  introduced.  For  single-user  or  frequency-division  multiplexed  wireless  systems,  we  de¬ 
velop  a  technique  we  refer  to  as  spread-response  precoding  which  replaces  the  interleaver 
traditionally  used  in  conjunction  with  coding  in  such  systems.  Prom  the  perspective  of 
the  coded  symbol  stream,  spread-response  precoding  effectively  transforms  an  arbitrary 
Rayleigh  fading  channel  into  a  nonfading,  simple  white  marginally  Gaussian  noise  chan¬ 
nel  with  no  intersymbol  interference.  Furthermore,  spread-response  precoding  requires  no 
additional  power  or  bandwidth,  and  is  attractive  in  terms  of  computational  complexity, 
robustness,  and  delay  considerations. 

In  the  multiuser  case,  spread-response  precoding  generalizes  to  a  new  class  of  orthogonal 
code-division  multiple-access  (CDMA)  systems  for  efficient  communication  in  environments 
subject  to  multipath  fading  phenomena.  The  key  characteristic  of  these  new  systems, 
which  we  refer  to  as  “spread-signature  CDMA”  systems,  is  that  the  associated  signature 
sequences  are  significantly  longer  than  the  interval  between  symbols.  Using  this  approach, 
the  transmission  of  each  symbol  of  each  user  is,  in  effect,  spread  over  a  wide  temporal  and 
spectral  extent,  which  is  efficiently  exploited  to  combat  the  effects  of  fading. 

Both  efficient  signature  sets  and  efficient  receiver  structures  for  such  systems  are  devel¬ 
oped.  Several  aspects  of  the  performance  of  the  resulting  spread-signature  CDMA  systems 
are  presented,  including  both  the  achievable  bit  error  rate  characteristics  as  well  as  the 
effective  capacity  of  such  systems.  The  results  suggest  that  spread-signature  CDMA  may 
be  an  attractive  alternative  to  conventional  CDMA  in  a  variety  of  application  scenarios. 
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Part  I 


Spread-Response  Precoding 


I.l  Introduction 

The  need  to  reliably  transmit  analog  and  digital  data  over  channels  subject  to  fading  arises 
in  a  wide  range  of  applications  including  mobile  radio  and  personal  wireless  systems,  and 
audio  and  television  broadcasting.  Generally,  the  fading  characteristics  of  the  channel, 
which  are  a  function  of  both  the  nature  of  the  transmission  media  and  the  relative  motions 
of  the  transmitter  and  receiver  in  a  given  system,  lead  to  variations  in  the  quality  of  channel 
both  in  time  and  in  frequency. 

Diversity  techniques  are  widely  used  in  communication  systems  to  compensate  for  these 
variations  [1]  [2].  These  range  from  simple  multiple  transmission  strategies  in  time,  fre¬ 
quency,  and  space,  to  more  sophisticated  diversity  techniques  based  on  the  use  of  coding. 
In  such  scenarios,  coding  is  used  to  combat  both  the  effects  of  fading  and  the  effects  of 
stationary  additive  noise. 

In  order  for  coding  to  be  effective  against  fading  in  particular,  it  is  generally  necessary 
to  combine  coding  with  interleaving,  a  simple  but  nevertheless  useful  form  of  precoding. 
The  purpose  of  interleaving  is  to  scramble  the  coded  data  stream  so  that  fading  channel 
effectively  seen  by  this  stream  is  uncorrelated  from  time-sample  to  time-sample.  This  sub¬ 
stantially  reduces  the  coding  complexity  required  to  achieve  a  given  level  of  fidelity,  allowing 
shorter  lengths  in  the  case  of  block  codes,  or  fewer  states  in  the  case  of  convolutional  codes. 

In  this  part  of  the  report,  we  develop  an  attractive  alternative  to  interleaving  which  we 
term  “spread-response  precoding.”  With  spread-response  precoding,  the  fading  chaimel  as 
seen  by  the  coded  data  stream  is  effectively  transformed  into  a  simple  additive  white  noise 
channel.  As  a  result,  when  combined  with  coding  techniques  such  as  trellis-coded  mod¬ 
ulation,  the  preceding  stage  combats  any  fading  effects,  while  the  coding  stage  combats 
the  remaining  additive  noise.  This  partitioning  appears  to  be  rather  attractive  in  terms 
of  system  complexity  considerations.  Indeed,  spread-response  precoding,  in  requiring  com¬ 
paratively  simple  linear  signal  processing  at  transmitter  and  receiver,  is  significantly  less 
demanding  computationally  than  error-correcting  coding  and  decoding  algorithms.  More¬ 
over,  precoding  constitutes  a  diversity  strategy  that  incurs  no  additional  cost  in  terms  of 
bandwidth  or  power,  and  is  competitive  with  traditional  approaches  in  terms  of  robustness 
and  delay  considerations. 

While  it  has  generally  been  understood  in  the  communications  community  that  inter- 
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x[n] 


x[n] 


Figure  1.1:  Precoding  system. 


leaving  is  not  generally  the  most  efficient  of  precoding  strategies  (see,  e.g.,  [3]),  the  literature 
has  offered  surprisingly  few  alternatives  to  date.  Perhaps  the  work  closest  in  spirit  to  the 
ideas  presented  here  is  contained  in  [4],  although  in  this  case  the  author  developed  a  form 
of  precoding  for  use  in  conjunction  with  interleaving  rather  than  as  an  alternative  to  inter¬ 
leaving. 

This  part  of  the  report  is  organized  as  follows.  In  Section  1.2,  we  outline  the  Rayleigh 
fading  channel  model  we  consider.  In  Section  1.3,  we  develop  the  transmitter  portion  of 
the  system,  and  discuss  fundamental  design  and  implementation  tradeoffs.  In  Section  1.4, 
we  turn  our  attention  to  the  receiver  portion  of  the  system,  developing  optimum  designs. 
In  Section  1.5,  we  explore  aspects  of  the  performance  of  systems  using  precoding,  including 
theoretical  capacity  and  attainable  bit  error  rate  calculations.  Finally,  Section  1.6  contains 
some  concluding  remarks. 

1.2  System  Model 

Fig.  I.l  depicts  a  block  diagram  of  the  overall  system  we  consider.  In  this  figure,  x[n]  is 
the  complex-valued  M-ary  symbol  sequence  representing  the  coded  bit  stream,  and  y[n] 
is  the  precoded  symbol  stream  to  be  transmitted.  The  transmitted  data  y[n]  is  corrupted 
by  complex-valued  fading  and  additive  noise,  producing  r[n]  at  the  receiver.  The  received 
data  is  first  processed  by  an  equalizer  to  produce  y[n],  then  by  a  postcoder  to  produce  x[n]. 
Finally,  a  decoder  (not  shown)  processes  :r[n]  to  produce  an  estimate  of  the  original  bit 
stream. 

The  channel  in  Fig.  1.1  is  the  equivalent  discrete-time  baseband  model  of  a  fairly  general 
stationary  Rayleigh  fading  channel  with  uncorrelated  scattering  and  bandwidth  Wo-  The 
channel  consists  of  two  components,  a  linear  time-varying  filter  which  captures  the  effects 
of  multipath  fading  due  to  multiple  scatters  in  the  transmission  medium,  and  an  additive 
noise  term  representing  both  receiver  noise  and,  more  significantly,  co-channel  interference. 
More  specifically,  the  response  of  the  channel  to  an  input  sequence  y[n]  is  given  by 

r[n]  =  5]  o[n;  k]  y[n  —  k]  -I-  i(;[n]  (I.l) 

k 

where  w[n]  is  a  zero-mean  complex  stationary  white  Gaussian  sequence  with  variance 

E  [l«;[n]p]  =  TVoWo, 

and  a[n;k],  the  response  of  the  channel  at  time  n  to  a  unit-sample  at  time  n  —  k,  is  a 
complex  Gaussian  fading  process.  For  fixed  values  of  k,  the  o[n;  k]  are  zero-mean  complex 
jointly  stationary  and  Gaussian  sequences.  Furthermore,  uncorrelated  scattering  implies 
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sequences  corresponding  to  distinct  values  of  k  are  statistically  independent.  Hence, 


E  [a[n;  k]  a*[n  —  m;  i]]  =  Ra[m;  k]  5[A;  -  1] 
where  5[n]  is  the  unit-sample,  ie., 

j-r  ,  A  1  1  n  =  0 

=  I  0  otherwise  ’  ('-2) 

With  uncorrelated  scattering,  the  time-variant  channel  frequency  response^ 

A{uj;  n]  =  ^  o[n;  k]  (1.3) 

k 

is  then  stationary  in  both  n  and  lo  and  satisfies 

E[A{u\n^  =  0  (I.4a) 

E\[A{u}-,n]f^  =  al.  (I.4b) 

Both  w[n]  and  a[n;  k]  are  assumed  to  be  statistically  independent  of  the  input  to  the  channel. 

Finally,  we  assume  that  while  the  transmitter  does  not  have  access  to  the  fading  channel 
kernel  a[n;  k]  or  its  statistics,  these  parameters  are  known,  or,  more  typically,  can  be  reliably 
measured  at  the  receiver. 


1.3  Spread- Response  Precoding 

In  this  section,  we  consider,  specifically,  the  transmitter  portion  of  the  system  in  Fig.  I.l.  We 
begin  by  observing  that  although  the  detailed  characteristics  of  a  fading  channel  fluctuate 
from  time  sample  to  time  sample,  the  performance  of  communication  systems  using  such 
channels  is  generally  dictated  by  the  average  characteristics  of  that  channel  over  time. 
Certainly  most  capacity  estimates  for  such  channels  involve  averaging  of  this  type;  see,  e.g., 
[5]  and  the  references  therein.  As  a  consequence,  in  principle  an  efficient  communication 
strategy  for  such  channels  would,  in  some  sense,  “spread”  the  transmission  of  each  symbol 
over  a  large  nimiber  of  time  samples. 

Conveniently,  spreading  of  this  type  can  be  achieved  through  simple  linear  time-invariant 
(LTI)  filtering  of  the  coded  symbol  stream — a  form  of  precoding.  Specifically,  denoting 
the  miit-sample  response  of  the  precoding  filter  by  h[n],  the  transmitted  sequence  is^  (c/. 
Fig.  1.1) 

y[n]  =  x[n]  *  h[n]  =  ^  x[k]  h[n  —  k].  (1.5) 

k 


^We  adopt  the  useful  convention  of  using  peirentheses  (•)  to  denote  continuous-valued  arguments  and 
brackets  [•]  to  denote  discrete- valued  arguments.  For  functions  of  two  arguments  where  the  first  is  continuous 
and  the  second  is  discrete  (as  in  the  case  of  time-variant  frequency  responses)  we  use  the  convenient  mixed 
notation  (•;  •].  The  notation  [•;  •)  is  used  in  a  similar  manner. 

^We  use  operator  *  to  denote  convolution,  and  the  superscript  *  to  denote  complex  conjugation. 
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For  convenience,  we  restrict  our  attention  to  the  case  in  which  /i[n]  is  a  real-valued  sequence. 

It  is  also  highly  desirable  for  the  precoding  to  be  lossless,  in  which  case  y[n]  constitutes 
an  orthonormal  transformation  of  the  data  symbols  a:[n].  Lossless  LTI  filters  satisfy  the 
time-domain  constraint 

h\n  —  k\  h[n  —  ^]  =  —  ^]-  (L6) 

n 

In  the  firequency  domain,  the  condition  (1.6)  corresponds  to 

=  1,  (1-7) 

and,  for  this  reason,  lossless  filters  are  frequently  referred  to  as  allpass  filters  in  the  signal 
processing  literature  [6].  An  important  property  of  such  filters  which  follows  immediately 
from  (1.6)  is  that 

=  h[— n]. 

Accordingly,  we  may  conveniently  and  stably  recompute  a;[n]  from  y[n]  according  to 

=  y[n]  *  h[-n]  =  y[^]  h[k  —  n].  (1.8) 

k 

From  a  practical  standpoint,  it  is  generally  necessary  to  restrict  our  attention  to  finite 
impulse  response  (FIR)  precoders.  However,  as  is  well  known,  the  only  lossless  FIR  filters 
are  the  shifted  unit-samples,  i.e.,  h[n]  =  <i[n  — fc]  for  arbitrary  k.  Nevertheless,  many  infinite 
impulse  response  (HR)  lossless  filters  have  strongly  localized  temporal  support  and  can  be 
therefore  be  truncated  without  significantly  altering  their  characteristics.  More  generally, 
there  are  a  wide  variety  of  FIR  filters  that  closely  approximate  the  losslessness  condition 
(1.6). 

To  make  this  notion  of  approximation  more  precise,  consider  the  class  of  LTI  filters  h[n] 
with  unit-energy,  i.e., 

=  1,  (1.9) 

n 

and  let 

=  /i[n]  *  h[—n]  =  h[k]  h[k  —  n] 
k 

denote  the  autocorrelation  function.  Then,  since  the  autocorrelation  of  a  perfectly  lossless 
filter  is  the  unit-sample,  a  useful  measme  of  the  deviation  from  losslessness  is  the  reciprocal 
of  the  total  sidelobe  energy  in  the  autocorrelation  function.  Specifically,  we  may  define  the 
following  convenient  losslessness  merit  factor  for  a  unit-energy  filter  fi[n]: 

C"  =  (e  (®oN  -  ■*["1)")  =  (2  E  w)  (I-IO) 

with  large  corresponding  to  filters  that  are  nearly  lossless.  More  specifically,  we  have 
0  <  <  00  with  the  right-hand  equality  if  and  only  if  the  filter  is  perfectly  lossless.  This 

merit  faetor  is,  in  fact,  identical  to  that  introduced  by  Golay  to  evaluate  the  quality  of  low 
autocorrelation  binary  sequences  [7],  and  we  note  that  it  is  sometimes  useful  to  exploit  an 
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equivalent  frequency-domain  definition  for  viz., 

which  can  be  verified  by  applying  Parseval’s  theorem  to  (I.IO). 

While  losslessness  is  an  important  attribute  of  a  precoding  filter,  a  second  important 
attribute  concerns  their  effectiveness  in  their  primary  function:  spreading  the  transmission 
of  each  symbol  over  a  large  number  of  time  samples.  This  is  achieved  when  the  precoder’s 
unit-sample  response  energy  is  widely  dispersed  in  time,  or,  equivalently,  when  the  precoding 
system  has  strong  partial  response  characteristics.  Indeed,  one  can  interpret  spread-response 
precoding  as  a  form  of  partial-response  precoding,  although  the  objectives  of  traditional 
partial-response  precoding  are  markedly  different. 

A  useful  measure  of  dispersion  for  an  arbitrary  unit-energy  filter  h[n]  in  the  context  of 
this  work  is  given  by 

with  large  corresponding  to  good  spreciding  characteristics.  Using  (1.9),  it  follows  that 
for  unit-energy  filters 

Vh  >  1  (1.12) 

with  equality  when  h[n\  =  <5[n].  However,  for  FIR  filters  of  length  N  with  unit-energy,  (1.9) 
also  implies  that 

Vh<N 

with  equality  precisely  when 

|h[n]|  =  l/\/jV,  all  0  <  n  <  JV  -  1. 

Consequently,  for  FIR  precoders  maximum  dispersion  {i.e.,  Vh  =  N)  is  obtained  when  h[n] 
is  an  antipodal  (binary)  sequence.  Conveniently,  binary  sequences  are  also  highly  attractive 
in  terms  of  both  computational  efficiency  and  numerical  sensitivity. 

While  both  good  losslessness  and  spreading  characteristics  are  desirable  in  design  of  FIR 
precoders,  they  are  competing  objectives.  Based  on  the  preceding  discussion,  at  one  extreme 
the  precoder  h[n]  =  (5[n]  corresponds  to  the  best  possible  but  the  worst  possible  Vh-  At 
the  opposite  extreme,  precoders  with  binary-valued  unit-sample  responses  provide  the  best 
possible  T>h  for  a  given  length  constraint,  but  poor  values  of  C^.  Indeed,  a  conjecture  of 
Golay  [8]  based  on  an  ergodicity  postulate  suggests  that  for  such  binary  sequences 

max£^  — >  12.3247 •  ■  •  as  N  oo. 

h 

Consequently,  as  we  will  discuss  later,  for  a  fixed  filter  length  prescribed  by  external  delay 
constraints  imposed  on  the  overall  system,  the  precoder  design  process  requires  a  compro¬ 
mise  between  the  two  components.^ 


®  Actually,  as  we’ll  discuss  in  Section  II.5.1,  this  problem  can  be  circumvented  by  using  only  a  slightly  more 
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When  spread-response  precoding  is  used  in  a  system,  the  transmitted  data  have  some 
rather  special  asymptotic  characteristics.  First,  we  note  that  with  lossless  precoding  we  get, 
via  (1.7),  that  the  transmitted  stream  y[n]  will  have  same  power  spectrum  as  the  original 
coded  data  x[n],  i.e., 

Syiu)  =  5x(a;). 

In  particular,  when  a:[n]  is  a  sequence  of  statistically  independent  complex-valued  symbols 
each  with  energy  Sg,  then  y[n]  is  a  complex  wide-sense  stationary  white  sequence  with  vari¬ 
ance  €s-  In  addition,  it  can  be  shown  that  in  the  limit  of  infinite  dispersion  (T>h  oo)  the 
transmitted  stream  is  marginally  Gaussian,  i.e.,  that  each  transmitted  sample  y[n]  has  a 
Gaussian  distribution.  This  follows  from  a  straightforward  Central  Limit  Theorem  argu¬ 
ment;  using  (1.5)  we  see  each  y[n]  is  the  balanced  sum  of  a  large  number  of  independent 
random  variables.  Prom  the  point  of  view  of  transmission  security  and  capacity  consider¬ 
ations,  such  characteristics  are  rather  appealing.  However,  we  should  also  note  that  from 
the  point  of  view  of  peak-to-average  power  and  receiver  synchronization  requirements,  these 
characteristics  present  important  practical  challenges  in  terms  of  system  design.  However, 
such  considerations  are  beyond  the  scope  of  the  present  report. 

1.3.1  Precoder  Design 

As  will  become  apparent,  the  overall  system  performance  depends  in  a  rather  complicated 
manner  on  the  loss  and  dispersion  characteristics  of  the  precoding  filter.  Furthermore,  even 
if  a  weighted  combination  of  the  losslessness  and  dispersion  factors  were  to  correspond  to 
some  reasonable  measure  of  ultimate  system  performance,  optimizing  such  a  criterion  is  an 
intrinsically  difficult  numerical  problem. 

Nevertheless,  there  are  a  number  of  approaches  that  work  weU  in  practice  for  obtaining 
precoding  filters  corresponding  to  a  reasonable  compromise  between  losslessness  and  dis¬ 
persion  factors,  and  providing  good  overall  system  performance.  In  fact,  design  techniques 
are  suggested  from  a  variety  of  sources,  largely  because  such  filters  find  application  in  a 
number  of  distinct  commimication  problenas  ranging  from  the  mitigation  of  impulse  noise 
in  communication  systems  [9]  to  robust  quantization  of  non-Gaussian  sources  [10].  Those 
we  discuss  here  are  closely  related  to  the  techniques  which  have  proven  useful  in  these 
problems. 

As  discussed  earlier,  binary  precoders  have  attractive  computational  properties  and  op¬ 
timal  dispersion  factors  but  poor  losslessness  factors.  The  binary  sequences  having  optimal 
losslessness  factors  for  a  given  length  constraint  N  are  tabulated'^  in  [11]  for  N  <  71.  For 
length  N  =  13,  the  corresponding  sequence  is,  for  example,  the  well-known  Barker  sequence 
for  which  14.08.  Likewise,  for  length  N  =  27,  the  optimum  sequence  has  w  9.85. 
Another  class  of  binary  sequences  with  better-than-average  losslessness  characteristics  are 
the  maximal  length  shift  register  sequences  (or  m-sequences)  [12],  for  which  ~  3  for 
large  N  as  shown  in  [13]. 


complicated  linear  precoder  structure — in  particular,  one  that  is  not  time-invariant  but  rather  periodically 
time- varying. 

^Actually,  to  simplify  the  search  for  such  sequences,  typically  only  skew-symmetric  binary  sequences 
are  considered.  This  restriction  is  not  severe  however,  since  this  class  generally  includes  many  of  the  best 
sequences  even  when  the  restriction  is  removed. 
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In  practice,  the  losslessness  characteristic^  of  these  sequences  can  be  markedly  improved 
while  largely  preserving  their  dispersion  proT:erties  by  relaxing  their  binary  amplitude  con¬ 
straints  in  a  controlled  maimer.  One  way  to  a  ’...omplish  this  is  to  apply  the  following  simple 
heuristic  algorithm.  Beginning  with  an  initial  binary  sequence  /loW?  we  proceed  to  maxi¬ 
mize  in  the  form  (I.IO)  via  an  iterative  ascent  algorithm  (such  as  the  Simplex  method) 
subject  to  the  unit-energy  constraint  (1.9).  Naturally,  a  global  optimization  would  result  in 
the  trivial  solution  <5[n].  However,  the  objective  function  typically  has  many  local  maxima, 
and  in  practice,  the  algorithm  generally  converges  to  such  a  fixed  point,  and  results  in  a 
filter  with  substantially  improved  at  a  relatively  modest  cost  in  V^. 

Other  filters  with  reasonable  dispersion  factors  are  also  suitable  as  an  initialization  for 
this  type  of  algorithm.  For  instance,  in  [10]  a  chirp  sequence  of  the  form 


fioM  =  /?sin 


7rn(n  —  1) 


forn  =  0, 1, . . . ,  —  1 


is  used  with  good  results.  This  initialization  is  especially  useful  in  the  design  of  very  long 
precoders,  in  which  case  the  corresponding  optimal  binary  sequences  are  unknown.  However, 
depending  on  the  specific  value  of  N  and  on  the  specific  iterative  descent  algorithm,  certain 
initializations  lead  to  better  precoders  than  others  in  terms  of  overall  system  performance. 
Consequently,  trial  and  error  is  invariably  involved. 


1.4  System  Characteristics  and  Receiver  Design 

In  this  section,  we  turn  our  attention  to  the  receiver  portion  of  the  system  in  Fig.  I.l.  We 
begin  by  noting  that  the  receiver  for  decoding  the  bit  stream  can  be  partitioned  into  two 
stages.  The  first  is  the  equalization  stage,  which,  as  depicted  in  Fig.  1.1  and  without  loss 
of  generality,  can  be  described  as  the  cascade  of  an  equalizer  and  postcoder.  The  second 
is  the  decoding  stage  (not  shown  in  Fig.  I.l),  typically  consists  of  some  form  of  Maximum 
Likelihood  (ML)  sequence  detection.  As  will  become  apparent,  when  no  coding  is  employed, 
simple  symbol-by-symbol  detection  generally  suffices  at  the  decoder.  However,  in  the  sequel 
we  restrict  our  attention  to  the  equalization  and  postcoding  stages. 

In  general,  the  equalizer,  which  compensates  for  the  fading,  is  a  linear  time- varying  filter 
whose  kernel  is  6[n;  A;],  so  that 


y\n\  =  ^  fe[n;  A:]  r[n  —  A:]. 
k 

Typically,  this  kernel  is  a  function  of  the  fading  channel  kernel  a[n;  A:]  and  the  noise  statistics, 
both  of  which  are  assumed  to  be  available  at  the  receiver. 

In  turn,  the  postcoder  inverts  the  transformation  of  input  symbols  that  takes  place 
during  precoding,  and  is  simply  a  linear  filter  whose  unit-sample  response  is  (c/.  (1.8))  a 
time-reversed  version  of  the  lossless  precoding  filter  h[n],  ie., 

x[n\  =  /i[— n]  *  y[n].  (L13) 

The  overall  system  consisting  of  the  channel  with  precoding,  equalization,  and  post- 
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coding,  we  refer  to,  for  convenience,  as  the  “composite  channel.”  In  the  remainder  of  this 
section  we  first  derive  the  key  properties  of  this  composite  channel,  then  optimize  them 
through  judicious  choice  of  the  equalizer  parameters. 

Our  main  result  is  that  subject  to  only  relatively  mild  ergodicity  constraints  on  channel 
the  use  of  lossless  precoding  with  a  large  dispersion  factor  leads  asymptotically  to  an  addi¬ 
tive  white  noise  composite  channel  that  is  not  only  free  of  fading,  but  has  no  inters3Tnbol 
interference.  In  order  to  make  our  result  precise,  we  first  define  a  sufficiently  realistic  class 
of  ergodic  kernels  for  our  purposes,  which,  for  convenience,  we  term  “admissibly  ergodic.” 

Definition  1  Let  /[n;  A:]  he  the  kernel  of  a  linear  system,  and  define 

f[n-,k]  =  f[n-,k]- E[f[n-,k]].  (1.14) 

Furthermore,  let 

d[n;  k]  =  Y^  f[n\  1]  f*[n  -hj-k] 
l 

and  define 

d[n-,  fe]  =  d[n;  k]—E  [d[n;  A:]] . 

Then  /[n;  A:]  is  an  admissibly  ergodic  kernel  if  the  following  conditions  are  satisfied: 


E[f[n-,k]]  =  /x5[A:]  for  every  k,n  (I.17a) 

E^f[n;k]f*[n  —  m-,l]^  =  J?[m;  A:]  5[A:  —  i]  for  every  k,l,m,n  (I.17b) 

E  j^d[n;  A:]  d*[n  —  m;/]j  =  T[m;k,l]  for  every  k,l,n,m  (I.17c) 

Sr  =  ^^|J?[n;A;]|  <  oo  (I.17d) 

k  n 

St  =  {I-17e) 

k,l  m 


Before  proceeding,  we  adopt  some  convenient  nomenclature.  In  general  it  will  be  conve¬ 
nient  to  view  a  generic  linear  kernel  such  as  /[n;  A:]  as  a  collection  of  sequences  in  n  indexed 
by  k.  Hence,  when  we  refer  to  “the  sequence  f[n;  k]”  we  specifically  mean  the  sequence  in 
n  corresponding  to  a  fixed  (but  generally  arbitrary)  value  of  k.  Prom  this  viewpoint,  the 
conditions  in  Definition  1  are  straightforward  to  interpret.  Conditions  (I.17a),  (I.17b),  and 
(1. 17c)  are  essentially  stationarity  constraints.  They  ensure,  specifically,  that  the  kernel  se¬ 
quences  /[n;  A:]  and  the  correlation  sequences  d[n;  A:]  are  each  jointly  wide-sense  stationary. 
The  condition  (1. 17b)  also  ensures  that  sequences  f[n;k]  corresponding  to  distinct  values 
of  k  are  imcorrelated.  Finally,  conditions  (I.17d)  and  (I.17e)  in  effect  ensure  that  linear 
combinations  of  the  kernel  sequences  f[n;  A;]  are  mean-  and  correlation-ergodic. 

Equivalently,  the  conditions  in  Definition  1  can  be  interpreted  in  terms  of  stationarity 
and  ergodicity  constraints  on  the  time-variant  system  frequency  response 

F(u;;n]  =  '^f[n;k]e-^‘^\  (1.18) 

k 

We  note,  in  particular,  that  (I.17a)  and  (I.17b)  imply  that  F(a;;n]  is  wide  sense-stationary 


(1.15) 

(1.16) 
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in  both  n  and  u  and  satisfies 


E  n]]  =  jjL. 


(1.19) 


Furthermore,  with 

F’(w;  n]  =  n]  —  E  [F’(w;  n]] 

we  have 

E  [f(0;  n]  F*(0  -  w;  n  -  m]]  =  $(u;;  m],  (1.20) 

where  ^’(a;;m]  is  the  spaced-frequency  spaced-time  correlation  function  of  the  system  and 
satisfies 

(w;  m]  =  ^2,  k]  (1-21) 

k 

Finally,  for  completeness,  we  note  that 

S(A;  fe]  =  53  (1-22) 

m 

is  the  system’s  scattering  function, 

r(A)  =  53  ^(0;  m]  (1.23) 

m 

is  the  system’s  Doppler  power  spectrum,  and 

cri  =  i?[0;  A:]  =  var  /[n;  k]  (1-24) 

is  the  multipath  intensity  profile  or  delay  power  spectrum  of  the  system.  Hence,  the  total 
power  is  given  by 

=  var  [F(w;  n]]  =  $(0;  0]  =  53  (1-25) 

k 

We  can  now  present  our  main  theorem  concerning  the  composite  system  depicted  in 
Fig.  I.l.  A  proof  is  presented  in  Appendix  III.A. 

Theorem  1  Let  a:[n]  be  a  sequence  of  zero-mean,  uncorrelated  symbols,  each  with  energy 
£s;  let  a[n;  A:]  and  u;[ra]  be  as  defined  in  (I.l);  and  let  c[n;A:]  denote  the  kernel  of  the  com¬ 
posite  linear  system  formed  by  cascading  the  channel  corresponding  to  kernel  a[n;  A:]  with 
the  equalizer  corresponding  to  kernel  &[n;A;],  i.e., 

c[n;  fc]  =  53  ~  (1.26) 

l 

Finally,  suppose  c[n;  A;]  and  h[n\  k]  are  both  admissibly  ergodic  kernels  in  the  sense  of  Defi¬ 
nition  1.  Then,  as  Vh  oo,  we  have^,  for  each  n, 

x[n]  -h  u[n],  (1-27) 

where  v\n\  is  a  complex-valued,  marginally  Gaussian,  zero-mean  white  noise  sequence,  un- 

®We  use  the  notation  to  denote,  specifically,  convergence  in  the  mean-square  sense. 


12 


correlated  with  the  input  symbol  sequence  rc[n]  and  having  variance 


varri[n]  =  Sscr^  +  AfoVVoial  +  iMbl^)- 


(1.28) 


Effectively,  Theorem  1  asserts  that  the  use  of  sufficient  precoding  transforms  the  channel 
“seen”  by  the  coded  symbol  stream  from  a  fading  channel  into  a  marginally  Gaussian 
white  noise  channel  with  no  intersymbol  interference.  As  a  result,  the  characteristics  of  the 
composite  channel  depend  only  on  the  statistics  of  the  fading  channel  parameters  a[n;  k] 
and  ^t;[n]  but  not  on  the  values  of  the  parameters  themselves.  Not  surprisingly,  this  is  a 
natural  consequence  of  the  time  averaging  induced  by  precoding  being  applied  to  an  ergodic 
fading  process. 

A  few  additional  remarks  regarding  Theorem  1  are  appropriate.  First,  as  can  be  readily 
verified  from  the  proof  in  Appendix  III.A,  Theorem  1  is  in  fact  true  even  when  the  kernel 
a[n;  k]  and  receiver  noise  u;[n]  are  not  Gaussian.  This  observation  is  important  in  terms  of 
robustness,  since  these  processes  are  of  course  at  most  approximately  Gaussian  in  practice. 

Second,  we  emphasize  that  v\p\  in  Theorem  1  is  a  marginally  Gaussian  process.  Specif¬ 
ically,  this  means  that  v[n]  is  a  sequence  of  uncorrelated  random  variables  each  having  a 
Gaussian  distribution,  but  that  the  random  variables  are  not  necessarily  jointly  Gaussian. 
Hence,  u[n]  is  not  necessarily  a  Gaussian  random  process,  and,  as  a  result,  there  may  exist 
at  least  some  statistical  dependence  among  the  noise  samples.  However,  we  shall  assume 
that,  in  practice,  u[n]  is  an  at  least  approximately  Gaussian  process,  so  that  the  statistical 
dependency  among  samples  is,  in  some  sense,  negligible. 

The  noise  v[n]  has  other  special  characteristics  as  well.  As  is  apparent  from  (1.28), 
the  noise  consists  of  the  sum  of  two  uncorrelated  components.  The  first  component  h£^ 
power  A/oWo(cr^  -t-  \pbf)  and  is  due  to  the  noise  in  the  original  fading  channel.  The  second 
component  has  power  and  is  inherently  generated  in  the  precoding  process.  The 
existence  of  this  second  noise  component  means  that  boosting  the  transmitter  power  in  the 
system  also  leads  to  an  increase  in  noise  power  in  the  composite  channel.  For  this  reason, 
we  shall  find  that  the  familiar  matched-filter  equalizer  [14]  is  not  best  suited  for  use  with 
precoding. 

Finally,  we  point  out  that  although  Theorem  1  estabhshes  only  an  asymptotic  result 
valid  for  perfectly  lossless  precoding  with  infinite  dispersion  (implying,  for  example,  infi¬ 
nite  delay),  the  asymptotic  behavior  can  be  approximated  arbitrarily  closely  with  realizable 
precoders.  In  particular,  with  suitably  chosen  finite  length  precoders,  the  white  marginally 
Gaussian  channel  model  is  an  excellent  approximation  to  the  actual  composite  channel. 
Furthermore,  the  delay  requirements  in  this  case  are  comparable  to  those  of  conventional 
interleaving.  Specifically,  the  filter  length  N  required  to  effectively  converge  to  the  equiva¬ 
lent  model  is  given  by 

iV  «  AT'  (Ta  +  1)  (1.29) 

where  N'  is  the  length  required  in  the  case  of  memoryless  fading,  and  where  is  the 
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coherence  time  of  the  fading  channel  (in  samples).  Specifically, 


1  00  /  1  °°  \ 

I]  ’^a(0;m]  =  -2  1]  -  1 

o  m=l  \  ®  m=0  ) 


(1.30) 


where,  defined  according  to  (1.20),  ^a(a;;m]  is  the  original  channel’s  spaced-frequency 
spaced-time  correlation  fimction.  We  can  alternatively  express  in  terms  of  the  (double¬ 
sided)  Doppler  spread  of  the  fading  channel  which  we  define  as 


/OO 

Ta(A), 

-OO 


Ta(0) 


(1.31) 


where,  defined  according  to  (1.23),  Ta(A)  is  the  Doppler  power  spectrum  of  the  channel.  In 
particular,  comparing  (1.31)  with  (1.30)  we  see  that  our  definitions  of  coherence  time  and 
Doppler  spread  are  related  according  to 

1  -f-  2Ta  =  ‘l-nj  be,. 

Later  we  discuss  appropriate  values  for  N'  in  practice. 

It  is  important  to  note  that  while  large  coherence  times  (or,  equivalently,  small  Doppler 
spreads)  correspond  to  larger  inherent  delays  in  the  system,  they  need  not  incur  addi¬ 
tional  computational  complexity  either  at  the  transmitter  or  receiver.  In  fact,  from  a  single 
prototype  precoder  h[n]  one  can  derive  an  entire  family  of  precoders  with  the  same  loss,  dis¬ 
persion,  energy,  and  computational  characteristics  but  suitable  for  scenarios  corresponding 
to  different  coherence- time/Doppler-spread  parameters.  Specifically,  if  fi[n]  is  a  precoder 
designed  for  memoryless  fading  (r^  =  0),  then  the  corresponding  precoder  firaW  for  the 
case  in  which  the  fading  has  coherence  time  is  obtained  by  simply  upsampling  fi[n],  ie.. 


h[n/M\  n  = 


-M,  0,  M,  2M, 


otherwise 


(1.32) 


where® 


M  =\Ta  +  1] 


is  the  upsampling  factor.  Equivalently,  this  upsampling  relationship  can  be  described  in 
the  frequency  domain  as 

=  H{Mu}).  (1.33) 

Using  both  (1.32)  and  (1.33)  one  can  verify  that  the  energy,  loss,  dispersion,  and  computa¬ 
tional  characteristics  are  unaffected  by  such  upsampling. 

Let  us  now  consider  the  design  of  the  equalizer.  We  begin  by  observing  that  Theorem  1 
implicitly  imposes  certain  constraints  on  the  equalizer  kernel  &[n;  fc]  in  order  that  the  equiv¬ 
alent  channel  model  structure  is  attained  asymptotically.  These  constraints  are,  from  a 
practical  standpoint,  in  fact  relatively  mild.  Specifically,  it  suffices  that  b[n;  fc]  be  chosen  so 


^The  ceiling  function  fx]  denotes  the  smallest  integer  greater  than  or  equal  to  x. 
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that  both  b[n;  fc]  and  the  cascade  of  a[n;  k]  with  6[n;  A:],  (ie.,  (1.26))  are  admissibly  ergodic. 
For  convenience,  let  us  refer  to  kernels  6[n;  k]  with  this  property  as  admissible  equalizers. 

Among  the  class  of  admissible  equalizers,  some  yield  better  composite  channels  than 
others.  One  useful  measure  of  the  quality  of  the  composite  channel  is  the  signal-to-noise  ratio 
(SNR).  Certainly  when  u[n]  is  Gaussian  the  both  the  theoretical  capacity  and  achievable  bit 
error  rates  increase  monotonically  with  the  channel  SNR.  Consequently,  a  useful  criterion 
for  equalizer  design  is  to  select  among  the  admissible  equalizers  that  yielding  the  largest 
SNR  in  the  composite  channel.  Conveniently,  when  6[n;  k]  is  an  admissible  equalizer,  the 
SNR  in  the  composite  channel  follows  directly  from  Theorem  1  as 


j(b)  = 


IMcI^ 

<^c  +Co(ff6  +  |M6P) 


(1.34) 


where 

Co=A’oWo/£*.  (1.35) 

In  Sections  1.4.1  and  1.4.2,  we  proceed  to  derive  these  optimum  equalizers  in  the  two  cases 
of  perhaps  greatest  interest  in  practice,  corresponding  to  frequency-nonselective  fading  and 
frequency-selective  slow  fading,  respectively. 


1.4.1  Prequency-Nonselective  Fading 

In  this  section,  we  restrict  our  attention  to  the  Rayleigh  fading  channel  model  of  (I.l)  with 

a[n;  k]  =  a[n]  5[A:]. 

in  which  case,  as  is  straightforward  to  show,  the  admissible  equalizers  are  also  of  the  form 

6[n;  k]  =  6[n]  5[A:]. 


Accordingly,  we  may  rewrite  (1.34)  in  this  scenario  as 


7(6)  = 


\E[ab]\‘^ 

var  [a6]  +  0^1^]  ’ 


(1.36) 


where  ^0  is  as  defined  in  (1.35)  and  where  we  have  omitted  a  specification  of  the  time  sample 
n  due  to  stationarity. 

To  derive  the  optimum  equalizer,  we  begin  by  rewriting  (1.36)  as 


7(6)  = 


l/(p{b)  -  1 


where 


£;[(!ap  +  eo)m’ 

and  where  we  have  exploited  the  identity 


(1.37) 


(1.38) 


var[a6]  =  E  pa6|^j  —  |F[a5]p 
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where  the  (complex)  constant  of  proportionality  is  arbitrary. 

Not  surprisingly,  the  optimum  equalizer  is  specified  only  to  within  an  arbitrary  gain 
factor,  since  such  factors  do  not  affect  the  resulting  SNR.  In  addition,  it  is  interesting  to 
note  (and  can  be  readily  verified)  that  the  equalizer  yielding  optimum  SNR  in  the  composite 
channel,  ie.,  (1-39),  also  corresponds,  when  suitably  normalized,  to  a  minimum  mean- 
square  error  linear  equalizer  for  the  fading  channel.  Specifically,  y[n]  and,  in  turn,  x[n]  are 
minimum  mean-square  error  linear  estimates  of  y\n]  and  x[n],  respectively.  Although  he 
did  not  establish  its  optimality,  it  was  Wittneben  [4]  in  his  preliminary  work  on  precoding 
for  the  nonselective  fading  channel  with  interleaving  who  first  suggested  that  a  minimum 
meain-square  error  type  equalizer  was  well-suited  to  this  scenario. 


1.4.2  Frequency-Selective  Fading 

In  this  section,  we  consider  the  more  general  channel  model  (LI)  for  which  the  fading  is,  in 
general,  frequency  selective.  This  scenario  generally  arises  when  larger  transmission  band- 
widths  are  used,  as  is  often  desirable  to  achieve  additional  diversity  benefit.  Furthermore, 
in  deriving  the  corresponding  optimum  equalizers,  we  exploit  the  fact  that  in  this  case  the 
fading  process  becomes  increasingly  slowly-varying  as  the  bandwidth  is  expanded. 

Accordingly,  we  begin  by  assuming  that  the  coherence  time  in  the  channel  is  large,  so 
that  the  fading  coefficients  are  effectively  constant  over  a  several  time  samples.  Then  the 
time-variant  channel  frequency  response  A(a;;n]  will  vary  slowly  with  n,  and  admissible 
equalizers  will  also  have  a  slowly- varying  time- variant  frequency  response  B{ui\  n].  Further¬ 
more,  we  have 

C(a;;n]  w  A(ci;;n]R(a;;n].  (1-40) 

Now  (1.19)  and  (1.25)  imply 

Pc  =  F[C(u;;n]] 

(Jc  =  var[C(a;;n]] 

and 

+  \y-b?  =  E  [|S(a;;  njl^]  .  (I.41c) 

Hence,  using  (1.41)  with  (1.40)  in  (1.34),  we  get  that  the  channel  SNR  is  effectively  given 


(I.41a) 

(I.41b) 
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by 


7(6) 


\E[AB]\ 


(1.42) 


var[AS]  +  ^o^[|5i2]’ 

where,  again,  ^0  is  as  given  by  (1.35),  and  where  we  have  omitted  specification  of  both  the 
time  sample  n  and  frequency  value  ui  due  to  stationaxity.  Clearly,  (1.42)  is  identical  in  form 
to  (1.36),  and,  consequently,  the  equalizer  providing  the  best  SNR  for  the  composite  channel 
follows  analogously  as 

A*{u};n] 


B{u;;n]  oc 


|A(a;;n]|2  +^o 


(1.43) 


where,  again,  the  (complex)  constant  of  proportionality  is  arbitrary.  As  in  the  nonselective 
fading  case,  we  remark  that  this  equalizer  not  only  maximizes  the  channel  SNR  but  also, 
when  suitably  normalized,  makes  x[n]  a  minimum  mean-square  error  linear  estimate  of  x[n]. 

The  receiver  structure  in  this  frequency-selective  scenario  warrants  some  additional  dis¬ 
cussion.  With  b[n;  fe]  varying  slowly  with  n,  let  us  denote  by  b[k]  the  nominal  value  of  b[n;  fc] 
over  some  sufficiently  long  generic  time  interval.  To  implement  the  unit-sample  response 
b[k],  it  is  useful,  at  least  conceptually,  to  view  the  equalization  process  in  two  stages.  The 
first  stage  implements  the  numerator  of  (1.43),  and  corresponds  to  the  appropriate  matched- 
filter  equalizer,  and,  hence,  is  a  conventional  RAKE  receiver  [1].  The  unit-sample  response 
of  this  stage  is,  therefore,  a*[— fc].  The  second  stage,  which  implements  the  denominator  of 
(1.43),  then  performs  additional  compensation,  taking  into  accoxmt  the  SNR  in  the  channel. 
The  unit-sample  response  corresponding  to  this  second  stage,  which  we  denote  by  e[k],  is 
symmetric  and,  in  general,  infinite  in  extent  (hence,  two-sided),  precluding  any  recursive 
implementation.  However,  typically,  the  tails  of  the  unit-sample  response  fall  off  quickly.  In 
particular,  when  o[A:]  is  non-zero  for  only  finitely  many  values  of  A:  as  is  firequently  assumed 
in  practice,  e[A:]  decays  exponentially  quickly  with  k.  Hence,  using  truncation  e[/!:]  may  be 
effectively  approximated  as  a  symmetric  FIR  filter.  As  an  example.  Fig.  1.2  depicts  the 
normalized  RMS  value  of  the  coefficient  e[fc],  viz., 


v'BBeli)P)/£)l|e[01P]. 

as  a  function  of  k,  where 


E 


1/2  k  =  0,l 
0  otherwise 


(1.44) 


and  where  ^0  =  a^. 


1.5  Performance 

In  this  section,  we  explore  the  potential  performance  achievable  through  the  use  of  spread- 
response  precoding  with  optimized  receivers.  We  begin  by  observing  that  with  the  optimum 
equalizer,  i.e.,  (1.39)  in  the  case  of  firequency  nonselective  fading  or  (1.43)  in  the  case  of 
frequency  selective  fading,  we  readily  obtain  that  the  corresponding  SNR  in  the  composite 
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Figure  1.2;  A  typical  RMS  unit-sample  response  of  the  second  stage  of  equalizer.  In  this 
example,  the  variances  of  the  fading  coefBcients  a[A:]  were  chosen  according  to  (L44),  and 
the  SNR  was  chosen  so  that  =  <^0- 
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channel  is  given  by 


where® 


|A(a;;n]p£:5 


(1.45) 


(1.46) 


AToWo  ^  ^  ^ 

denotes  the  SNR  in  the  original  fading  channel  at  a  particular  time  instant  n  and  frequency 
w.  A  useful  notion  of  the  capacity  of  the  composite  channel  is  given  by  the  equivalent 
Gaussian  capacity^,  i.e.,  using  (1.45), 


C/Wo  =  -iog(l  +  7o)  =  -log  (e 

\  L70  +  1J/ 


Exploiting  the  readily  verified  identity 


where 


E  [log(7o  +  1)]  =  [:^]  =  £;x(Co) 


(1.47) 


(1.48) 


(1.49) 


denotes  the  average  SNR  in  the  original  fading  channel  (c/.  (I.4b)),  and  £^i(-)  denotes  the 
exponential  integral  [15] 


l(^)  =  f 

Ju 


00  g-t 


we  get  that  (1.51)  can  be  expressed  more  conveniently  as 


C/Wo  =  -log(Coe<‘>Ei(Co)). 


(1.50) 


(1.51) 


The  capacity  estimate  (1.51)  can  be  compared  to  some  related  capacity  calculations.  In 
particular,  by  the  Schwarz  inequality  we  get 

E  1  =  1 

[70  +  1]  E  [70  +  1]  1  +  1/Co  ’ 

so  that,  as  expected,  (1.51)  is  upper  bounded  by  the  capacity  of  the  Gaussian  channel  or, 
equivalently,  the  fading  channel  with  infinite  spatial  diversity,  ie.. 


Coo/Wo  =  log  (1  +  E  [70])  =  log  (1  +  1/Co)  • 


(1.52) 


More  generally,  there  have  been  a  variety  of  attempts  to  estimate  the  capacity  of  fading 
channels  in  the  literature  [5].  As  an  example,  between  C  and  Coo  lies  the  estimate  of  fading 


®Of  course,  in  the  case  of  nonselective  fading  A(u;;n]  =  a[n;  0]  =  a[n]  using  the  notation  of  Section  1.4.1. 
^More  precisely,  this  can  be  interpreted  as  a  constrained  capacity — specificaJly  the  bit  rate  that  can  be 
achieved  when  the  remauning  higher-order  statistical  dependencies  in  the  composite  channel  model  are  not 
exploited. 
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SNR  (dB) 


Figure  1.3:  Capacity  estimates,  Rayleigh  fading  channel.  The  solid  curve  is  the  capacity 
estimate  C  determined  using  precoding.  For  comparison,  the  dotted  curve  is  the  capacity  Coo 
when  infinite  spatial  diversity  is  available,  while  the  dashed  curve  is  the  capacity  estimate 
Co  due  to  Lee  [16]. 

channel  capacity  without  spatial  diversity  derived  by  Lee  [16],  i.e., 

Co/Wo  =  E  [log  (1  +  7o)]  =  ^i(Co),  (1-53) 

where  the  second  equality  follows  from  (1.48).  In  Fig.  1.3,  these  capacity  estimates  (in 
bits/sec/Hz)  are  plotted  as  a  function  of  the  average  available  SNR.  The  capacities  C,  Co, 
and  Coo  are  represented  by  the  solid,  dashed,  and  dotted  curves,  respectively. 

At  high  SNR  (i.e.,  small  Co),  we  can  use  the  series  expansion  [15] 


where  To  =  0.57721  •  •  •  is  Euler’s  constant,  to  show  that  (1.51)  is  given  asymptotically  by 


~  log 


log(l/Co)> 


By  contrast,  in  the  same  regime  (1.52)  behaves  like 

Coc/Wo  ~  log(l/Co) 

while  (1.53),  as  can  be  shown  via  (1.54),  satisfies 

Cq/Wq  ~  To  +  log(l/Co)- 


(1.55) 


(1.56) 


(1.57) 


Hence,  comparing  (1.57)  with  (1.56)  and  (1.55)  we  can  verify  that  at  high  SNR  the  difference 
between  Cqo  and  Cq  is  approximately  2.506  dB  while  the  difference  between  Cq  and  C  diverges 
to  infinity. 

We  stress  that  these  capacity  estimates  correspond  to  the  case  in  which  the  transmitter 
has  no  knowledge  of  the  state  of  the  faxling  channel  or  its  statistics  at  any  point  in  time, 
i.e.,  there  is  no  side-channel  for  feedback  from  the  receiver  to  the  transmitter.  We  note, 
however,  because  of  memory  in  the  fading  channel,  the  availabihty  of  a  feedback  path  would 
naturally  lead  to  higher  capacity  [17]. 

We  also  remark  that  to  approach  the  capacity  C  of  the  composite  charmel  requires,  of 
course,  that  coding  be  applied  to  the  data  stream  prior  to  precoding.  However,  because 
the  composite  channel  is  effectively  an  additive  white  Gaussian  noise  charmel,  any  of  the 
traditional  forms  of  coding  for  this  channel  would  be  appropriate.  In  particular,  we  note 
that  conventional  implementations  of  trelhs-coded  modulation  appear  to  be  well-suited  to 
this  scenario. 

Nevertheless,  even  without  coding  significant  improvements  in  bit  error  rate  performance 
can  be  achieved  from  the  inherent  diversity  benefit  of  using  spread-response  precoding. 
We  emphasize  that  this  is  in  marked  contrast  to  the  use  of  interleaving,  which  offers  no 
improvement  in  bit  error  rate  performance  without  coding.  For  the  purposes  of  illustration, 
let  us  consider  the  case  in  which  x[n]  is  an  uncoded  QPSK  (quadrature  phase-shift  keying) 
stream.  When  precoding  is  used,  the  bit  error  probabihty  as  a  function  of  the  SNR  per  bit. 


is  given  by 


where 


Mom  2Mom' 


(1.58) 


1 


and  where  7,  via  (1.45)  with  (1.46),  is  given  by 


Coei'^EiiCo)  ^ 
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with  (o  given  by  (1.49).  For  comparison,  without  precoding  the  QPSK  bit  error  probability 
using  L-fold  spatial  diversity  and  maximal  ratio  combining  is  given  by  [1] 


ro[L]  =  -  1 


Furthermore,  (1.59)  specializes  to 


[-TOS©(dfci)‘l  ™ 


when  there  is  no  spatial  diversity,  and  to 


nioo]  =  Q  (v^)  =  Q  (l/^/^) 


(1.60) 


when  there  is  infinite  spatial  diversity.  In  Fig.  1.4,  we  plot  bit  error  probability  as  a  function 
of  SNR  per  bit  (i.e.,  l/(2<^o))  with  and  without  precoding.  The  solid  curve  in  Fig.  1.4 
corresponds  to  the  use  of  precoding  (with  infinite  dispersion)  but  no  spatial  diversity.  The 
dashed  curves  correspond  to  the  use  of  no  precoding  but  T-fold  diversity  for  £  =  1,  2,  . . . ,  5. 
Finally,  the  dotted  curve  corresponds  to  the  use  of  no  precoding  but  infinite  spatial  diversity, 
Zj  — ^  oo. 

Comparing  V  with  'Po[l]  we  see  that  precoding  markedly  improves  bit  error  rate  per¬ 
formance  in  the  channel.  We  can  further  show,  by  applying  (1.54)  and  the  asymptotic 
expansion  [15] 


(—!)"» (2m)!  ^_2m 
y/^v  ^  2”*  ml 


(1.61) 


to  (1.58),  that  at  high  SNR  {i.e.,  small  Co)  the  bit  error  rate  with  precoding  is  given  by 


/log(l/Co) 


l/(2Co) 


By  contrast,  in  the  high  SNR  regime,  Vq[L]  is  well-approximated  as  [1] 


(1.62) 


^o[£]  ~  Co"", 


while  (1.60),  via  (1.61),  takes  the  form 

Po[oo]  ~ 

Thus,  we  note  that  the  while  V  falls  off  at  a  slower  rate  than  'Pofoo],  it  falls  off  faster 
than  ■Po[£]  for  any  fixed  L.  This  implies  that,  asymptotically,  precoding  provides  higher 
“effective  diversity”  than  can  be  achieved  using  spatial  diversity  with  any  L  <  oo,  but  less 
than  can  be  achieved  with  infinite  spatial  diversity. 

While  the  solid  curve  in  Fig.  1.4  is  obviously  a  lower  bound  on  the  bit  error  rate  per¬ 
formance  achievable  in  practice  using  finite  length  precoders,  it  is  quite  realistic  when  the 
precoding  filter  has  length  N  given  by  (1.29)  for  N'  ~  100.  In  Fig.  1.5,  we  plot  bit  error 


SNR  (dB/bit) 


Figure  1.4;  Bit  error  probabilities  using  uncoded  QPSK  on  the  Rayleigh  fading  channel. 
The  solid  curve  represents  the  performance  achievable  with  precoding.  For  comparison, 
also  depicted  is  the  attainable  performance  without  precoding  but,  instead,  with  spatial 
diversity  and  maximal  ratio  combining.  The  dotted  curve  corresponds  to  the  performance 
with  infinite  spatial  diversity,  while  the  successively  lower  dashed  curves  correspond  to  the 
performance  with  L  =  1,  2,  . . . ,  5  branches  of  spatial  diversity,  respectively. 
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probability  for  uncoded  QPSK  as  a  function  of  SNR  per  bit  using  realizable  precoders. 
The  solid  curve  indicates  the  infinite  dispersion  bound,  while  the  successively  lower  broken 
curves  indicate  the  performance  with  practical  precoders  corresponding  to  N'  =  27  and 
N'  =  200,  respectively. 

Finally,  we  remark  that  experiments  exploring  the  sensitivity  of  system  performance  to 
the  parameters  of  the  equalizer,  while  beyond  the  scope  of  this  report,  are  clearly  warranted 
in  the  future.  In  the  meantime,  preliminary  simulations  reported  in  [4]  involving  minimum 
mean-square  error  type  equalizers  lend  at  least  some  insight  into  what  kind  of  behavior 
may  be  expected.  Specifically,  the  results  suggest  that  performance  ought  to  depend  only 
weakly  on  the  estimate  of  Sg/  (A/o^o)  used  in  the  equalizer,  and  that  the  performance  can  be 
expected  to  degrade  fairly  gracefully  with  errors  in  estimates  of  the  magnitude  and  phase 
of  the  fading  coefficients. 


1.6  Summary 

Spread-response  precoding  as  developed  in  this  part  of  the  report  constitutes  a  potentially 
attractive  alternative  to  interleaving  in  a  wide  range  of  communication  systems  designed  for 
use  with  multipath  fading  channels.  Even  when  no  additional  coding  is  used,  precoding  can 
significantly  improve  system  performance  over  other  uncoded  systems.  Similarly,  the  use 
of  precoding  in  conjimction  with  coding  has  the  potential  to  substantially  reduce  computa¬ 
tional  complexity  at  both  the  transmitter  and  the  receiver  for  a  given  level  of  performance. 
This  is  because  the  effects  of  fading  are  entirely  controlled  by  the  precoding,  which  requires 
only  low-complexity  signal  processing.  Thus,  only  additive  noise  remains  for  coding  to  con¬ 
trol,  which  can  in  turn  be  achieved  with  comparatively  shorter  codes.  As  an  additional 
potential  featme,  the  noise-like  characteristics  of  the  transmitted  stream  resulting  from  the 
use  of  precoding  appears  to  be  well-suited  for  applications  involving  secure  conomunication. 

Nevertheless,  several  technical  issues  remain  to  be  explored.  For  example,  a  detailed 
investigation  of  the  complexity  benefits  that  can  be  realized  through  the  use  of  precoding 
with  coding  is  clearly  warranted.  The  development  of  techniques  for  adequately  synchro¬ 
nizing  the  receiver  in  such  systems  as  well  as  for  managing  the  inherent  peak-to-average 
power  requirements  is  also  an  important  avenue  for  future  research.  Furthermore,  in  order 
for  precoding  to  be  effective,  reliable  estimates  of  the  fading  channel  parameters  and  statis¬ 
tics  must  be  available  at  the  receiver.  It  appears  that  this  can  be  accomplished  in  practice 
through  the  use  of  training  sequences  or  decision  feedback  techniques.  However,  such  issues 
remain  to  be  explored  in  detail.  Finally,  as  remarked  earlier,  regardless  of  the  technique 
used,  any  such  parameter  estimates  are  imperfect  and,  as  a  result,  sensitivity  analysis  would 
ultimately  be  important  in  further  verifying  the  viability  of  spread-response  precoding. 
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Bit  Error  Probability  Pr(e) 


Figure  1.5:  Bit  error  probabilities  using  uncoded  QPSK  on  the  Rayleigh  fading  channel 
with  realizable  precoders.  The  dash-dotted  and  dashed  curves  represent  the  performance 
obtained  using  FIR  precoders  with  N'  =  27  and  N'  =  200,  respectively.  The  solid  curve 
represents  the  ideal  precoder  performance  bound  (N'  ooj. 
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Part  II 


Spread-Signature  Code-Division 
Multiple  Access 

II.  1  Introduction 

Systems  for  efficiently  coordinating  communication  among  multiple  users  in  multipath  fad¬ 
ing  environments  are  important  in  a  wide  range  of  applications.  Indeed,  such  systems  are 
essential  in  proposed  digital  mobile  radio  communications,  personal  wireless  systems,  indoor 
wireless  networks,  and  digital  audio  and  television  broadcasting  systems.  However,  rapidly 
escalating  demand  for  both  wider  availability  of  such  services  and  increasingly  sophisticated 
capabilities  has  put  great  pressure  on  the  limited  available  bandwidth  within  the  radio  spec¬ 
trum.  Given  such  constraints,  it  is  clear  that  the  use  of  increasingly  sophisticated  signal 
processing  in  wireless  modems  will  be  critical  to  accommodating  large  numbers  of  services 
and  users  within  the  available  spectrum. 

In  Part  I,  we  explored  the  use  of  spread-response  precoding  for  mitigating  the  effects  of 
fading  in  single-user  or  frequency-division  multiplexed  data  transmission  systems  used  in 
multipath  fading  environments.  Prom  the  perspective  of  the  coded  symbol  stream,  this  pre¬ 
coding  effectively  transforms  a  fairly  general  Rayleigh  fading  channel  into  a  nonfading,  sim¬ 
ple  white  marginally  Gaussian  additive  noise  channel  with  no  intersymbol  interference.  By 
using  such  precoding  to  combat  fading  while  reserving  coding  to  combat  only  the  remaining 
additive  noise,  substantial  reductions  in  system  complexity  appear  possible.  Furthermore, 
although  spread-response  precoding  represents  a  form  of  time-diversity,  it  is  efficient  in  the 
sense  that  it  requires  no  additional  power  or  bandwidth. 

In  this  part  of  the  report,  we  develop  a  natural  generalization  of  the  precoding  con¬ 
cept  for  general  multiuser  communication  problems  in  multipath  fading  environments.  The 
result  is  a  code-division  multiple-access  (CDMA)  system  in  which,  in  effect,  precoding  is  em¬ 
bedded  directly  into  each  user’s  signature  sequence  while  maintaining  orthogonality  among 
users.  We  term  the  resulting  system  “spread-signature  CDMA.”  These  signature  sequences 
have  some  very  special  characteristics  as  we  will  show,  most  notable  of  which  is  that  their 
temporal  extent  significantly  exceeds  the  intersymbol  (baud)  duration. 

In  a  manner  analogous  to  spread-response  precoding,  using  such  signature  sets  in  mul¬ 
tipath  fading  environments  has  the  effect  of  transforming  the  collection  of  channels  seen  by 
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the  individual  symbol  streams  from  a  collection  of  coupled  Rayleigh  fading  channels  into 
a  uncorrelated  collection  of  identical  nonfading  simple  white  maxginally-Gaussian  additive 
noise  channels.  In  essence,  spread-signature  CDMA  converts  various  degradations  due  to 
fading,  co-channel  interference,  and  receiver  noise  into  a  single,  comparatively  more  benign 
form  of  uncorrelated  additive  noise  that  is  white  and  quasi-Gaussian. 

Transformations  of  this  type  are,  in  general,  highly  desirable  in  multiuser  systems;  see, 
e.g.,  [18].  As  we  will  see,  spread-signature  CDMA  systems  provides  some  important  poten¬ 
tial  performance  advantages  over  traditional  CDMA  systems.  Furthermore,  we  will  see  that 
spread-signature  CDMA  systems  achieve  this  benefit  without  requiring  additional  power  or 
bandwidth,  and  are  attractive  in  terms  of  computational  complexity,  robustness,  and  delay 
considerations. 

Our  model  consists  of  a  fairly  general  cellular  multiple-access  scenario  in  which  each 
cell  contaiiK  a  single  base  station  (sometimes  referred  to  as  a  “cell  site”)  and  a  number 
of  mobiles  (or,  more  generally,  “subscribers”).  We  assume  that  both  forward  link  (base- 
to-mobile)  and  reverse  link  (mobile-to-base)  communication  is  required,  but  takes  place  on 
separate  (ie.,  non-interfering)  channels.  Between  each  transmitter-receiver  pair  is  a  fairly 
general  Rayleigh  fading  channel,  which  may  be  frequency  selective  or  non-selective. 

The  outline  of  the  Part  II  is  as  follows.  In  Section  II.2,  we  develop  the  equivalent  discrete¬ 
time  baseband  model  for  our  multiuser  system  and  some  basic  notation.  In  Section  II.3, 
we  then  develop  a  useful  framework  for  characterizing  the  generalized  orthogonal  CDMA 
signature  sets  of  interest  in  this  work,  and  use  this  framework  to  develop  an  efficient  family 
of  spread-signature  sets.  Section  II.4  develops  some  important  aspects  of  spread-signature 
CDMA  systems,  ranging  from  transmission  characteristics  to  channel-transformation  prop¬ 
erties  of  such  systems.  We  then  exploit  these  properties  to  develop  efficient  spread-signature 
CDMA  receivers.  Section  II.5  explores  the  potential  performance  of  such  optimized  systems 
both  in  terms  of  capacity  and  bit  error  rate  characteristics,  and  Section  II.6  contains  some 
concluding  remarks. 

II. 2  System  Model 

Consider  a  single  cell  of  a  multiple-access  system,  in  which  there  are  M  users,  all  sharing 
a  total  fixed  bandwidth  MWo,  so  that  Wo  is  the  efiective  bandwidth  per  user.  In  the 
equivalent  discrete-time  baseband  model  for  the  system,  the  modulation  process  can  be 
viewed  as  follows.  The  coded  symbol  stream  of  the  mth  user  {l<m<M),  which  we  denote 
by  Xm\p\i  is  modulated  onto  a  unique  signatme  sequence  hm[ri\  to  produce  ym\p\  which  is 
transmitted  within  the  total  available  bandwidth. 

Conceptually,  it  is  convenient  to  view  the  modulation  process  in  two  stages.  As  depicted 
in  Fig.  II.l,  these  stages  correspond  to  upsampling  (i.e.,  zero-insertion)  by  a  factor  M, 
followed  by  linear  time-invariant  filtering  with  the  signature  sequence,  i.e., 

ym[n]  =  51  ~ 

k 

The  multiuser  channel  we  consider,  which  is  depicted  in  Fig.  II.2,  is  a  rather  general  sta¬ 
tionary  Raleigh  fading  environment  with  uncorrelated  scattering.  More  specifically,  am[n;  k] 
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Xm[n] 


ymfn] 


Figure  II.l:  Modulation  of  the  the  mth  user’s  coded  symbol  stream  Xm{n\  onto  a  signature 
sequence  hm[n]  for  transmission. 


represents  the  zeromean,  complex- valued  Gaussian  kernel  of  the  fading  channel  seen  by  that 
user.  Hence,  the  sequence  obtained  at  the  receiver  is 

®m[wj  J/m[^  d"  (II. 2) 

m  k 


where  w[n\  is  a  zero-mean,  complex- valued  stationary  white  Gaussian  sequence  with  variance^ 


E 


M^lP]  =A^o>Vo 


(II.3) 


that  is  statistically  independent  of  both  the  am[n;  A:]  and  the  ym  W-  In  general,  the  randomly 
time- varying  kernels  am[n;k]  capture  the  effects  of  multipath  fading  due  to  both  fluctua¬ 
tions  in  the  media  and  the  relative  motions  of  transmitters  and  receivers  in  the  system. 
Meanwhile,  w[n[  captures  both  receiver  noise  and  any  somrces  of  co-channel  interference  not 
otherwise  taken  into  account.  We  use  to  denote  the  (stationary)  time-variant  fre¬ 

quency  response  corresponding  to  the  mth  channel.  Our  notation  for  characterizing  other 
aspects  of  this  class  of  linear  randomly  time-varying  systems,  which  is  essentially  the  same 
as  was  used  in  Part  I,  is  summarized  for  convenience  in  Appendix  III.B. 

Two  special  cases  of  this  general  channel  model  are  of  primary  interest.  The  first, 
corresponding  to  base-to-mobile  transmission,  is  referred  to  as  the  forward  Unk.  In  this 
case,  the  messages  to  the  individual  mobiles  are  multiplexed  together  before  being  broadcast 
over  the  channel.  Prom  the  perspective  of  a  particular  receiver  all  messages  are  transmitted 
through  the  same  channel,  i.e.. 


ai[n;  k]  =  a2[n’,  k]  =  •  •  •  =  aAf[n;  k]  =  o[n;  k].  (II-4) 

The  second  case,  corresponding  to  mobile-to-base  transmission,  is  referred  to  as  the 
reverse  hnk.  In  this  case,  the  messages  are  transmitted  through  separate  channels  to  the 
base  station.  With  reasonable  physical  separation  among  mobiles,  as  we  will  assume,  the 
kernels  of  the  individual  mobile-to-base  channels  may  be  modeled  as  mutually  independent. 
We  shall  further  assume  that  no  base-to-mobile  feedback  channel  is  available  to  provide 
synchronization  information  to  the  mobiles.  Accordingly,  the  kernels  In  our  model 

capture  the  effects  of  both  multipath  fading  and  asynchronism  among  users. 


^Note  that,  for  convenience,  in  our  equivalent  model  the  channel  parameters  are  bandwidth-normalized, 
i.e.,  the  statistics  of  both  am[u;  A;]  and  w[n]  are  independent  of  the  bandwidth  parameter  M.  To  maintain 
the  proper  dependence  of  signal-to-noise  ratio  on  M,  we  therefore  bcindwidth-normalize  the  transmitted 
power  eis  well.  Indeed,  as  we  will  establish  in  Section  II.4.1,  when  the  coded  stream  a:m[n]  has  power  £m, 
the  transmitted  stream  ym[n]  has  power  Sm/M. 
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Figiire  II.2;  General  multiuser  fading  channel  model,  where  am[n‘,k]  denotes  the  randomly 
time-varying  linear  kernel  corresponding  to  the  mth  user. 

Finally,  we  make  the  reasonable  assumption  that,  in  both  forward  and  reverse  link 
cases,  while  the  transmitters  have  no  knowledge  of  the  channel  kernels  am[n',k]  or  their 
statistics,  these  parameters  are  known — or,  more  typically,  can  be  reliably  measured — at 
the  corresponding  receivers  in  the  system. 

II.3  Orthogonal  Multiuser  Modulation 

In  traditional  CDMA  systems,  the  signature  sequences  hm[n]  used  in  the  modulation  (II.  1) 
have  length  equal  to  the  upsampling  rate  or  intersymbol  period  M.  In  this  way  the  signa¬ 
tures  are  used  in  a  nonoverlapping  manner.  However,  in  this  section  we  consider  signatures 
of  arbitrary  length  N.  When  N  >  M,  the  resulting  sequences  have  valuable  partial  response 
characteristics  akin  to  those  of  spread-response  precoding,  and  we  refer  to  these  sequences 
as  “spread-signature”  sequences. 

A  useful  mathematical  framework  for  representing  such  signature  sets  arises  out  of  mul¬ 
tirate  system  theory  [19],  as  we  now  show.  To  begin,  we  first  express  the  signature  set  as 
the  vector  sequence, 

-  h[n]  =  [ /ii[n]  h2[n]  •••  j  ,  (II.5) 

and,  for  convenience,  let  us  restrict  our  attention  to  real-valued  signature  sets.  When  each 


^The  superscript  ^  denotes  transposition. 
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of  the  component  signatures  has  only  finitely  many  non-zero  values,  we  shall  refer  to 
the  signature  set  as  having  finite-spread.  Specifically,  when 

h[n]  =  0,  n  <  0,  n  > 

we  say  that  the  signature  set  has  temporal  spread  N. 

Although  the  theory  can  accommodate  more  general  classes  of  signature  sets,  we  restrict 
our  attention  to  those  satisfying  certain  convenient  orthogonality  conditions,  which  facilitate 
both  analysis  and  implementation.  Specifically,  we  require  that  the  signature  sequences 
together  with  all  translates  by  integer  multiple  of  M  constitute  an  orthonormal  set,  ie.,^ 

h[A:  —  nM]  —  mM]  =  J[n  —  m]  I  (11.6) 

k 

where  I  denotes  the  identity  matrix  of  appropriate  size,  and  where  <5[n]  denotes  the  unit 
sample,  viz., 

t-r  1  A  I  1  n  =  0  /TT 

“  I  0  otherwise  ' 

In  turn,  the  corresponding  completeness  condition  for  this  orthonormal  set  can  be  expressed 
as 

h^[n  —  kM]  h[m  —  kM]  =  ^  hi[n  —  kM]  hi[m  —  kM]  =  5[n  —  m],  (II.8) 

k  k,i 

and,  in  fact,  (II.8)  can  be  interpreted  as  a  special  instance  of  Mercer’s  Theorem  [20]. 

It  is  important  to  emphasize  that  orthonormal  signature  sets,*.e.,  sets  satisfying  (II.6) 
and  (II.8),  correspond  to  lossless  systems.  As  a  result,  demodulation  in  the  absense  of 
distortion  is  particularly  simple.  To  illustrate,  if  the  input  symbol  streams  are  and 

the  superimposed  output  is  the  sequence  y[n],  i.e.,  via  (II.  1), 

2/N  =  H  ymW  =  hm[n  -  Mk]  Xm[k],  (II.9) 

^  m,k 

then  the  sequences  Xm[n]  can  be  reconstructed  from  y[n]  via 

Xm[n]  =  Y  ~  y  W-  (II.IO) 

k 

As  we  shall  see,  the  discrete-time  matched-filter  and  downsample  operation  (II.IO)  is  also 
a  key  component  of  the  demodulation  process  when  distortion  is  present  as  weU. 

For  M  >  2,  there  is  a  rich  collection  of  signature  sets  satisfying  (II.6)  and  (II.8),  even 
when  we  restrict  our  attention  to  signatures  with  finite-spread.  For  M  =  1,  however,  the 
modulation  process  (II.  1)  is  equivalent  to  prefiltering  with  a  linear  time-invariant  filter  whose 
unit-sample  response  is  the  signature  sequence  hi[n].  In  this  case,  which  was  extensively 
developed  in  Part  I,  the  condition  (II.6)  is  equivalent  to  requiring  that  /ii[n]  be  an  allpass 
filter,  and  it  is  weU-known  that  non-trivial  finite-length  allpass  filters  do  not  exist. 


®Note  that  (II.6)  directly  incorporates  the  natural  requirement  that  each  signature  sequence  hra[n]  have 
unit  energy. 
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Several  aspects  of  orthogonal  signature  sets  are  more  conveniently  viewed  in  the  fre¬ 
quency  domain.  Accordingly,  we  express  the  set  of  Fourier  transforms  corresponding  to 
(II.5)  in  the  form 

/°®  •  A  r 

h[n]e-^^dw^[Hi{u)  ffsM  •••  FmH  ]  •  (II.ll) 

For  example,  the  orthogonality  properties  (II.6)  and  (II.8)  are  frequently  more  conveniently 
viewed  in  the  frequency  domain — as  an  illustration,  (IL8)  is  equivalent  to'* 

HtM  HM  =  ^  =  M,  (11.12) 

k 

which  is  sometimes  referred  to  as  a  (power)  complementarity  condition. 

More  importantly,  the  frequency  domain  representation  (II.ll)  leads  to  the  efficient 
factorization 

H(w)  =  Q(Ma;)A(a;),  (11.13) 

where  Q(a;)  is  a  square  matrix  and  A(w)  is  the  Fourier  transform  of  the  delay  chain  of 
order  M,  i.e., 

<5[n]  =  [  5[n]  5[n  — 1]  •••  5[n  -  M -h  1]  ]^ ,  (11-14) 

whence 

A(u;)  =  [l  •••  ]^.  (11.15) 

The  decomposition  (11.13)  is  referred  to  as  the  polyphase  representation  of  the  set,  and  Q(w) 
is  termed  the  associated  pol3q)hase  matrix.  As  is  true  for  multirate  systems  in  general  [19], 
polyphase  representations  of  signature  sequences  are  not  only  conceptually  useful,  but  lead 
to  computationally  efficient  modem  implementations  as  well. 

For  specifically  orthonormal  signatmre  sets,  the  associated  polyphase  matrix  satisfies  the 
special  property 

Q(£^)Qt(a;)=I,  (11.16) 

and  we  remark  that  matrices  satisfying  (11.16)  are  termed  paraunitary.  Furthermore,  it  is 
straightforward  to  verify  that  Q(w)  is  independent  of  u  if  and  only  if  the  signature  set  is 
not  spread  (i.e.,  N  =  M). 

The  polyphase  matrices  associated  with  some  familiar  orthogonal  signature  sets  provide 
useful  insight.  For  example,  the  polyphase  matrix  corresponding  to  time-division  multiple- 
access  (TDMA)  systems  is 

Q(a;)=I, 

while  that  corresponding  to  ideal  frequency-division  multiple-access  (FDMA)  systems  has 
(fc,  Z)th  element® 

0  <  w  <  TT. 

In  contrast,  for  discrete  Fomrier  transform  (DFT)  multiplexing,  Q(a;)  is  the  DFT  matrix, 

'‘The  superscript  *  denotes  the  conjugate-transpose  operation. 

®Note  that  Q(w)  is  both  conjugate  synmietric,  i.e.,  Q*(w)  =  Q(— w),  and  27r-periodic. 
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%•  6. 

[QHlw  = 

For  Hadamaxd  sequence  based  CDMA  systems,  for  which  we  will  develop  a  powerful 
generalization  in  Section  IL3.1,  we  have 


QM  =  s, 


where  S  is  the  Hadamard  matrix  of  appropriate  dimension.  Recall  that  the  Hadamaxd 
matrix  of  dimension  M,  viz.,  'Em,  where  M  is  a  power  of  two,  is  defined  recursively:  for 
M  =  2,  4,..., 

where  Si  =  1. 

CDMA  system  designers  are  frequently  interested  in  the  auto-  and  cross-correlation 
characteristics  of  the  signature  sequences.  Indeed,  the  autocorrelation  characteristics  gen¬ 
erally  affect,  for  example,  the  ability  of  a  receiver  to  synchronize  to  the  transmission,  while 
the  cross-correlation  characteristics  generally  affect  the  degree  and  nature  of  co-channel  in¬ 
terference.  Ideally,  therefore,  one  would  like  the  autocorrelation  of  each  signature  hk{inb[  to 
satisfy 

hk[n]  *  hfe[— n]  »  (5[n],  (11.17) 

and  the  cross-correlation  between  distinct  signatures  hk[n]  and  hi[ri\  to  satisfy 

hk[n]  *  hi[—n]  « 0  k^l.  (11.18) 

It  is  well-known  that  (11.17)  and  (11.18)  are  conflicting  objectives  for  traditional  signature 
sets  (see,  e.g.,  Welch  [21]  or  Saxwate  and  Pursley  [12]).  In  fact,  this  is  also  true  for  spread- 
signatmre  sets.  To  see  this,  let  us  define  quadratic  auto-  or  cross-correlation  merit  factors 
which  penalize  deviations  from  (11.17)  and  (11.18),  respectively.  Specifically,  analogous  to 
the  merit  factors  defined  by  Golay  [7],  let 

l/4i  =  ^  Hn^)\  -  <5[A:  -  if  div.  (11.19) 

Then,  using  (11.12),  it  is  straightforward  to  verify  that,  for  any  M  >  2, 


“M/2  “M/2 

Sm/2  -'^M/2 


— y"— 

tr  4/ 
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(M-1) 
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k  k,  Ij^k 

(M-1) 


du) 


du>  = 


M 


(11.20) 


Hence,  from  (11.20)  we  see  that  good  autocorrelation  characteristics  can  only  be  obtained  at 
the  expense  of  cross  correlation  characteristics,  and  vice  versa.  At  one  extreme  the  trivial 
signature  set  corresponding  to  TDMA  systems  has  perfect  autocorrelation  characteristics. 
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but  the  worst  possible  cross-correlation  characteristics®,  j.e., 


fh  _  J  oo  k  I 

I  1  ' 


At  the  opposite  extreme  we  have  the  signature  set  corresponding  to  ideal  FDMA  systems. 
This  set  has  perfect  cross-correlation  characteristics  but  poor  autocorrelation  characteris¬ 
tics;  in  particular, 


1/(M-1)  k  =  l 
oo  k^l  ' 


In  practical  CDMA  systems,  a  compromise  between  these  extremes  is  generally  sought. 

The  auto-  and  cross-correlation  characteristics  constitute  only  one  of  the  important 
issues  in  the  design  of  good  spread  signature  sets.  For  example,  it  is  also  important  that 
the  signature  sets  we  develop  are  effective  in  spreading  the  transmission  of  each  symbol 
of  a  user’s  transmission  over  a  large  range  of  time  samples  in  order  to  mitigate  the  effects 
of  fading.  This  is,  of  comse,  analogous  to  the  objectives  of  spread-response  precoding 
as  developed  in  Part  1.  Accordingly,  we  define  a  dispersion  factor  T>h  which  measures  a 
signature  set’s  spreading  capability  via 


1 


(11.21) 


where  V}^^  represents  the  dispersion  in  the  sequence  *-e., 

=  .  (11.22) 

Note  fi-om  (11.21)  that,  as  we  would  expect  from  any  reasonable  definition  of  dispersion,  the 
set  has  perfect  spreading,  i.e.,  D/j  — )■  oo,  if  and  only  if  every  signature  in  the  set  is  perfectly 
spread,  i.e.,  — >  oo  for  every  1  <  m  <  M. 

Important  insights  are  obtained  by  examining  what  values  D/j  can  take.  It  is  straight¬ 
forward  to  verify,  for  instance,  that  for  all  orthonormal  signature  sets 


Vh  >  1,  (11.23) 

and  that  this  bound  is  attained  when  h[n]  is  the  TDMA  signature  set.  More  importantly, 
at  another  extreme  we  have,  for  finite-spread  signature  sets  with  temporal  spread  N, 

Vh  <  N,  (11.24) 

with  equality  precisely  when,  for  each  1  <  m  <  M, 

l^mMI  =  l/v'^j  0  <  n  <  iV  —  1. 


®Iii  fa.ct,  as  can  be  verified  from  (11.19)  and  (11.20),  this  is  true  for  any  orthogonal  signature  set  whose 
sequences  have  allpass  Fourier  transforms. 
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Hence,  for  finite-spread  signature  sets,  maximum  dispersion  is  achieved  when  the  signature 
sequences  are  antipodal  (binary- valued).  For  this  reason,  we  refer  to  the  corresponding 
signature  sets  as  “maximally-spread.”  Note,  too,  that  because  they  are  discrete-valued, 
maximally-spread  signature  sets  are  especially  attractive  in  terms  of  computational  effi¬ 
ciency  and  numerical  stability. 

As  another  important  design  issue,  it  will  also  be  important  that  the  spread-signature 
sets  we  use  possess  what  we  refer  to  as  a  good  “partitioning”  characteristics.  In  particular,  as 
we  will  see,  good  partitioning  results  in  more  uniform  distribution  of  co-channel  interference 
among  users  in  the  system.  To  develop  this  concept,  we  define  the  following  modified 
correlation  function 

0A .  [n,  m]  =  53  [«  -  Mk]  hi  [m  -  Mk]  (11.25) 

k 

which  corresponds  to  correlating  hi\n]  with  a  version  of  fii[n]  in  which  all  but  every  Mth 
sample  is  replaced  with  zero.  While  the  complementarity  condition  (II.  8)  directly  implies 
that 

^hi [n,  m]  =  6[n  -  m],  (11.26) 

i 

this  condition  says  nothing  about  the  properties  of  each  of  the  M  terms  in  the  summation 
of  (11.26).  However,  as  will  become  apparent  in  Section  II.4,  the  set  has  good  partitioning 
characteristics  when  the  unit-sample  in  (11.26)  is,  in  some  sense,  distributed  uniformly 
among  the  M  modified  correlation  fimctions. 

To  malce  the  notion  of  good  partitioning  characteristics  more  precise,  we  let 

Qhi  [n,  m]  =  Qhi  rn]  -  (11.27) 

denote  the  deviation  from  ideal  partitioning  in  each  component.  It  is  straightforward  to 
show  that  Qhi  [^i  is  the  following  symmetry,  periodicity,  and  finite  energy  characteristics: 

©ft .  [n,m]=  Qhi  [m,  n] ,  (II.28a) 


0ft^.  [n,  m]  =  0ft.  [n  +  kM,  m  +  kM],  any  k, 


(II.28b) 


1/M2<53©2jn,m]<(l-1/M)2. 


(II.28c) 


A  signature  set  has  asymptotically  perfect  partioning  if  0ftjn,m]  can  be  made  arbitrarily 
small  using  sufficiently  long  signatures.  Accordingly,  we  define  the  following  partitioning 
factor: 


where 


Xh  M  ^  Xhi 


=  sup  53  ©I  Jn,m], 
Xhi  m  ^ 


(11.30) 


and  note  that  good  partitioning  corresponds  to  large  partitioning  factors,  and  vice  versa. 
It  is  important  to  note,  however,  that  signature  sets  with  good  dispersion  factors  do  not 
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necessarily  have  good  partitioning  characteristics^.  Nevertheless,  we  emphasize  that  we 
shall  be  primarily  interested  in  signature  sets  which  have  both  good  dispersion  and  good 
partitioning  characteristics. 

II.3.1  An  Optimum  Class  of  Spread-Signature  Sets 

In  this  section,  we  develop  a  family  of  orthogonal  signature  sets  that  are  optimal  in  the  sense 
of  being  maximally-spread,  i.e.,  having  the  best  dispersion  characteristics  for  a  given  length 
(or,  equivalently,  delay)  constraint.  As  we  discussed  earlier,  the  maximally-spread  condition 
is  achieved  precisely  when  the  signatme  sequences  are  binary-valued.  The  signatures  we 
construct  are  conveniently  obtained  out  of  a  powerful  paraunitary  generalization  of  the 
Hadamard  matrix,  and  have  a  computationally  efficient  recursive  synthesis  that  is  attractive 
in  terms  of  modem  implementation. 

Our  construction  is  based  on  the  polyphase  decomposition  of  a  signature  set  (11.13).  In 
particxilar,  we  rephrase  the  problem  of  designing  a  suitable  orthogonal  signature  set  H(w) 
into  problem  of  designing  a  smtable  paraunitary  polyphase  matrix  Q(a;).  Thus,  requiring 
that  the  desired  H(a;)  correspond  to  a  maximally-spread  signature  set  is  equivalent  to 
requiring  that  polyphase  matrix  sequence  q[n]  whose  transform  is  Q(tj)  be  binary-valued. 

We  begin  by  observing  that  the  Hadamard  matrix  S  is  one  matrix  satisfying  these 
properties.  Accordingly,  we  let  our  zero-th  order  polyphase  matrix  be® 

=  S,  (11.31) 


and  note  that  the  spread  of  the  corresponding  signature  set  is  N  =  M. 

To  obtain  signature  sets  for  which  N  >  M,  we  exploit  a  recursion  that  preserves  the 
binary  sequence  requirement;  specifically,  we  let 

Q«(a;)  =  2  A(M*-^u;)  Q^*-^)(u;),  i  =  1,  2,  . . .  (11.32) 

where  A(w)  is  the  diagonal  delay  matrix  whose  diagonal  is  constructed  from  the  elements 
of  A(u;),  i.e.,  with  A(a;)  as  defined  in  (11.15), 

A(u;)  =  diag  A(a;). 

We  see  immediately  that  the  paraunitary  property  (11.16)  is  preserved  by  the  recursion 
(11.32)  because  the  product  of  paraunitary  matrices  is  also  paraunitary.  It  is  similarly 
straightforward  to  verify  that  the  recursion  (11.32)  preserves  binary  sequence  property — 


^To  verify  this,  it  suffices  to  note  that  from  any  orthogonal  signature  set  we  can  construct  a  new 

orthogonal  signature  set  pm[n]  via 

9m[n]  =  ^  em[k]  hm[n  -  kM] 

k 

where  the  em[n]  are  any  set  of  nontrivial  lossless  (allpass)  filters.  However,  while  this  new  system  generally 
has  greater  dispersion  Vg  >  Vh,  it  is  straightforward  to  show  that  the  partitioning  factor  is  unchemged,  i.e., 
Xg  =  Xh- 

®For  convenience,  we  restrict  our  attention  to  orders  M  for  which  Hadamard  matrices  exist.  These 
include,  for  example,  all  integers  M  that  are  powers  of  two. 
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indeed,  in  the  time  domain  one  can  interpret  (11.32)  and  (11.31)  as  implementing  a  succession 
of  simple  sequence  concatenations  initiated  with  Hadamard  sequences. 

Using  (11.13),  we  can  also  express  the  recursion  (11.32)  directly  in  terms  of  the  signature 
set  vector  Fourier  transform;  specifically,  we  have,  for  i  =  1,  2,  . . 

HW(a;)  =  EA{M^oj)  H(*-^)(u;).  (11.33) 

Prom  (11.33)  we  can  verify  that  the  spread  of  the  signature  set  grows  by  a  factor  of  M  with 
each  application  of  the  recursion,  so  that,  in  particular,  has  spread  N  =  for  i  = 
0,  1,  2,  ....  For  convenience,  several  sets  of  signature  sequences  obtained  by  the  recursion 
(11.32)  with  (11.31),  and  corresponding  to  different  values  of  M  and  N,  are  tabulated  in 
Appendix  III.C. 

As  a  historical  aside,  it  is  interesting  to  note  that  orthogonal  systems  of  the  type  con¬ 
structed  in  this  section  have  been  reinvented  numerous  times  over  the  last  several  decades, 
particularly  for  the  case  corresponding  to  M  =  2.  However,  it  would  appear  that  this  work 
is  the  first  to  attempt  to  exploit  such  systems  in  multiuser  communication  problems. 

Using  a  variety  of  constructions,  such  systems  emerged  independently  in  a  variety  of 
unrelated  communities  within  mathematics,  physics,  and  engineering.  Work  within  the 
engineering  community  dates  traced  back  to  1950  when  Golay  constructed  pairs  of  pseu¬ 
dorandom  sequences  which  he  referred  to  as  “complementary  sequences”  [22]  [23].  These 
sequences — now  frequently  referred  to  as  Golay  sequences — were  defined  as  binary-valued 
sequences  satisfying  the  frequency-domain  complementarity  condition  (11.12)  for  M  =  2.^ 
These  Golay  sequences  have  subsequently  been  explored  extensively,  though  the  focus  has 
been  primarily  on  issues  of  existence  of  such  pairs  for  various  values  of  N;  see,  e.g.,  [24]  and 
[25]. 

Prom  (11.12)  it  is  apparent  that  binary-valued  orthogonal  sequences  are  a  subset  of  com¬ 
plementary  sequences.  However,  more  importantly,  complementary  sequences  also  useful  in 
the  construction  of  orthogonal  sequences,  which  both  Golay  and  Turyn  [26],  and  later  Taki, 
et  al.  [27]  observed  in  the  case^°  M  =  2.  The  corresponding  generalizations  for  M  >  2 
appear  in,  e.g.,  Tseng  and  Liu  [28]. 

Within  the  mathematics  community,  binary- valued,  pseudorandom,  orthogonal  sequence 
pairs  (and,  in  particular,  a  time-domain  version  of  the  recursion  (11.33)  for  M  =  2)  were 
also  discovered  independently  by  both  Shapiro  [29]  and,  later,  Rudin[30].  As  a  result,  the 
frequency  domain  representations  for  such  sequences  are  sometimes  referred  to  as  Rudin- 
Shapiro  polynomials. 

As  discussed  in,  e.g.,  Odlyzko  [31],  several  useful  properties  of  such  sequences,  as  well  as 
their  connections  to  other  families  of  binary  pseudorandom  sequences,  have  been  developed. 


®  As  such,  there  is  a  natural  correspondence  between  such  complementary  pairs  and  power  complementary 
filters  as  developed  in  the  raultirate  signal  processing  literature  [19]. 

^®To  see  this,  it  suffices  to  verify  using  (11.16)  with  (11.12)  that 


Q(w)  = 


A(ui)  B{u}) 
B*(w)  -A*(w) 


(11.34) 


is  a  paraunitary  matrix  whenever  the  corresponding  sequences  a[n]  and  6[n]  constitute  a  complementary 
pair. 
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For  example,  the  asymptotic  auto-  and  cross-correlation  characteristics  of  our  maximally- 
spread  orthogonal  signature  sets  for  the  case  M  =  2  follow  immediately  from  the  results  of 
Newman  and  Byrnes  [32],  In  particular,  we  have  that 

“^3,  N  CO,  (11.35) 

which,  using  (11.20),  implies  that  for  I, 

Ch  3/2,  N  CO.  (11.36) 

Based  on  earlier  discussion  in  Section  II.3,  (11.35)  and  (11.36)  suggest  that  our  maximally- 
spread  signature  sets  are  localized  in  neither  time  nor  frequency.  In  fact,  not  only  is  this  the 
case,  but  such  strong  spreading  in  both  time  and  frequency  is  critical  to  good  performance  in 
our  intended  application.  Furthermore,  in  addition  to  their  good  spreading  characteristics, 
our  maximally-spread  signature  sets  have  good  asymptotic  partitioning  characteristics  as 
well;  specifically,  Xh  oo  as  N  oo. 

II.4  System  Characteristics 

In  this  section,  we  develop  the  key  characteristics  of  spread-signature  CDMA  systems,  and 
then  apply  these  results  in  developing  efficient  receivers  for  such  systems.  We  begin  with  a 
general  discussion  of  transmission  characteristics. 

II.4.1  Transmission  Characteristics 

Spread-signature  CDMA  systems  give  rise  to  transmissions  with  some  rather  special  spectral 
and  temporal  characteristics,  only  some  of  which  are  shared  by  conventional  CDMA  systems. 
To  illustrate  the  spectral  features,  we  begin  by  observing  from  (II.  1)  that,  provided  the 
coded  symbol  stream  of  the  mth  user  is  stationary,  the  corresponding  transmitted  stream 
is  cyclostationary  with  a  time-averaged  power  spectrum  given  by 

where  Sx„{u})  is  the  power  spectrum  of  the  coded  symbols  of  the  mth  user.  When 
consists  of  uncorrelated  symbols  of  energy  Sm,,BS  we  shall  generally  assume  in  practice,  the 
transmitted  power  spectrum  associated  with  the  mth  user  further  simplifies  to 

SyM  =  \Hm{cu)\‘^£m/M.  (11.37) 

From  (11.37)  several  observations  can  be  made.  First,  we  immediately  note  that  since 
Hm{u!)  has  unit-energy  the  total  transmitted  power  in  this  case  is  Sm/M.  When  compared 
with  (II.3)  we  see  that  this  appropriately  results  in  a  signal-to-noise  ratio  (SNR)  that 
decreases  inversely  with  system  bandwidth. 

We  also  note  from  (11.37)  that,  in  general,  the  transmitted  power  spectrum,  is  not, 
strictly-speaking,  white.  However,  when  the  maximally-spread  signatme  sets  of  Section  II.3.1 
are  used,  the  transmitted  power  spectrum  is  broadband  and,  in  a  particular  sense,  asymp¬ 
totically  white.  More  specifically,  as  N  oo,  each  “converges”  to  a  highly  irregular 
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function  whose  energy  is  effectively  uniformly-distributed  over  frequency.  Furthermore,  the 
transmitted  spectrum  corresponding  to  forward  link  transmission  is  typically  truly  white. 
Indeed,  when  the  component  coded  symbol  streams  are  mutually  uncorrelated  and  of  equal 
energy  £m  =  £,  it  is  straightforward  to  verify  using  (11-12)  that  the  aggregate  transmitted 
power  spectrxun  is 

Sy{u})  =  £,  (11.38) 

where  y\n\  is  the  aggregate  transmission,  ie., 

yN  =  S]ymM-  (n.39) 

m 

The  temporal  characteristics  of  spread-signature  CDMA  transmissions  are  noteworthy 
as  well.  In  traditional  CDMA  systems  the  individual  (baseband)  transmissions  are  gener¬ 
ally  binary- valued.  However,  in  spread-signature  CDMA  systems,  although  the  signature 
sequences  and  coded  streams  may  be  binary-valued  (if,  e.g.,  the  signature  sequences  of  Sec¬ 
tion  II.3.1  are  used),  the  transmitted  streams  generated  via  the  modulation  (II.l)  generally 
are  not.  This  is,  of  comrse,  a  consequence  of  the  overlap  between  modulated  symbols  than 
invariably  results  from  choosing  signature  sequence  lengths  N  that  exceed  the  intersymbol 
epoch  M.  In  fact,  as  A"  oo,  a  simple  Central  Limit  Theorem  argument  suggests  that 
each  ym\p\  is  a  marginally  Gaussian  process. 

This  quasi- Gaussian  behavior  is  generally  rather  appealing  from  the  point  of  view  of 
certain  transmission  security  and  capacity  considerations.  However,  it  is  important  to  point 
out  that,  as  was  true  in  the  case  of  spread-response  precoding  systems,  these  characteristics 
also  pose  significant  engineering  challenges  in  terms  of  managing  peak-to-average  power 
and  receiver  synchronization  requirements.  While  such  issues  certainly  warrant  further 
investigation,  they  are  beyond  the  scope  of  the  present  report. 

11.4.2  Receiver  Characteristics  and  Design 

In  Section  II.4.1  we  explored  the  transmission  characteristics  of  spread-signatmre  CDMA 
by  examining  what  the  coded  symbol  streams  look  like  (after  modulation)  from  the  per¬ 
spective  of  the  channel.  In  this  section,  we  develop  efficient  receivers  for  such  systems  by 
examining  what  the  channel  looks  like  from  the  perspective  of  the  coded  symbol  stream 
(before  modulation). 

We  begin  by  observing  that  what  the  coded  symbol  stream,  in  effect,  “sees”  is  the 
cumulative  effect  of  modulation,  followed  by  distortion  introduced  by  the  channel,  followed 
by  processing  performed  at  the  receiver.  Accordingly,  we  refer  to  what  the  coded  symbol 
stream  sees,  then,  as  the  “composite  channel.” 

Before  we  explore  the  key  properties  of  this  composite  channel,  we  must  first  develop  an 
appropriate  receiver  structure  for  these  systems.  In  general,  the  receiver  processing  required 
to  recover  the  mth  transmitted  message  is  comprised  of  three  stages.  First,  the  received 
data  r[n]  is  equalized  according  to 

ym[n]  =  6m[n;  r[n  -  k].  (11.40) 

k 
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r[n] 


Figure  II.3:  Receiver  structure  for  extracting  the  symbol  stream  of  the  mth  user.  The  Srst 
stage  is  equalization,  producing  ym[n],  while  the  second  stage  is  demodulation,  producing 
XmM-  -A  hnal  stage  (not  shown)  is  decoding. 

As  is  apparent  from  (11.40)  and  as  is  depicted  in  Fig.  II.3,  the  corresponding  equalizer  is,  in 
general,  a  linear  time- varying  filter  whose  kernel  we  denote  by  bm[n;  A:].  Because  the  fading 
channel  coefficients  and  statistics  are  assumed  to  be  available  to  the  receiver,  the  equalizer 
kernel  6m[-;‘]  is  generally  a  function  of  all  the  channel  kernels  afc[-;-],  k  =  1,2, ...M, 
although  the  dependence  is  often  simpler  in  certain  special  cases  of  interest. 

In  the  second  stage,  also  depicted  in  Fig.  II.3,  the  equalized  data  is  demodulated  from 
the  corresponding  signature  sequence,  viz., 

~  y  ym  W  hm[k  Mil]. 

k 

This  is  conveniently  interpreted  as  a  discrete-time  matched-filter  and  downsample  operation. 

Finally,  the  last  stage  of  the  message  recovery,  which  is  not  depicted  in  Fig.  II.3,  con¬ 
sists  of  decoding  the  demodulated  stream  XmM  (using,  for  example,  maximiun  likelihood 
methods).  It  is  important  to  point  out,  however,  that  while  we  will  not  consider  the  actual 
implementation  of  the  decoder  in  this  report,  we  will  discuss  in  Section  II.5  what  kinds  of 
coding  and  decoding  strategies  are  appropriate  for  spread-signature  CDMA,  and  their  effect 
on  system  performance. 

The  composite  system  consisting  of  modulation,  the  channel,  equalization,  and  demod¬ 
ulation  has  some  appealing  characteristics  provided  the  channel,  and,  in  turn,  the  corre¬ 
sponding  equalizers,  have  some  reasonable  ergodicity  properties.  In  effect,  these  ergodic 
properties  are  required  to  ensure  that  spreading  the  transmission  of  eaeh  symbol  over  a 
sufficiently  long  time  interval  by  modulation  allows  these  symbols  to  see  the  average  char¬ 
acteristics  of  the  fading  channel — i.e.,  that  the  time-averaging  implicit  in  the  modulation 
is  equivalent  to  ensemble-averaging.  The  following  technical  definition  will  prove  to  be 
sufficient  for  om:  purposes. 

Definition  2  Let  /j[n;  k]  be  the  kernels  of  a  family  of  linear  systems,  and  define 

fj  [n;  fc]  =  fj  [n;  k]-E  [fj  [n;  A:]] .  (11-41) 

Furthermore,  let 

d)j  [n,  m;  k,  Z]  =  fj[n;  A:]  fp  [m;  Z]  (11.42) 

and  define 

d^j  [n,  m;  k,  Z]  =  [n,  m;  k,l]  —  E  ^dj  [n,  m;  k,  Z]  j  .  (11.43) 
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Then  the  family  fj[n;  A:]  is  admissibly  ergodic  if  the  following  conditions  are  satisfied: 


E 


E[fj[n-M 

E  [fj[n;k]fp[m-,l]j 
[n,  m;  k,  Z]  dj  [n' ,  m';  k\  Z']j 


=  E[Fj]S[k]  (IL44a) 

=  [n  —  m;  Ar]  5[A:  —  1]  (II.44b) 

=  Tp  [n  —  —  m',  n  —  n';  k,  I,  k',  Z']  (II.44c) 


n  k 

E  E  Tj  [ni,n2,n3;ki,k2,k3,ki] 


ki,k2 
"3  fc3,*4 


<  OO 


(IL44d) 

(II.44e) 


Definition  2,  in  fact,  represents  a  straightforward  generalization  of  the  corresponding 
definitions  constructed  in  Part  I.  In  particular,  conditions  (II.44a),  {IL44b)  and  (II.44c)  axe 
essentially  stationarity  constraints,  while  (II.44d)  and  (IL44e)  are  ergodicity  constraints. 
Note,  too,  that  for  convenience  we  have  omitted  specification  of  the  frequency  u;  and  time 
n  from  £^[fj(a;;n]]  in  (IL44a)  due  to  stationarity. 

The  key  properties  of  the  composite  channels,  which  we  now  develop,  represent  a  nat¬ 
ural  generalization  of  those  developed  for  single-user  systems  in  Part  I.  In  particular,  we 
show  that  subject  to  only  relatively  mild  ergodicity  constraints,  the  use  of  spread-signature 
modulation  with  sufficiently  long  signatures  leads  to  composite  channels  that  are  effectively 
a  collection  of  uncorrelated  additive  white  noise  channels,  each  of  which  is  not  only  free 
of  fading,  but  has  no  intersymbol  interference.  In  developing  our  detailed  results,  we  con¬ 
sider  the  forward  link  and  reverse  link  systems  separately  since  there  are  some  fundamental 
differences  in  behavior. 

The  Forward  Link 

The  characteristics  of  the  composite  system  for  the  forward  link  are  summarized  in  the 
following  theorem.  A  proof  is  provided  in  Appendix  III.D. 

Theorem  2  (Forward  Link)  Letxm[n\  be  mutually  uncorrelated  sequences  of  zero-mean, 
uncorrelated  symbols,  each  with  energy  6m,  o.nd  let  am[n;k]  and  t/;[n]  be  as  defined  in  (II. 2) 
and  (11-4) ■  Furthermore,  let  the  common  equalizer  be 


A 


Z)i[n;  A:]  =  b2\n]  A:]  =  •  -  •  =  A:]  =  fe[n;  A:], 


and  let  c[n;  A:]  denote  the  kernel  of  the  linear  system  formed  by  cascading  the  channel  whose 
kernel  is  a[n;  A:],  with  the  equalizer  whose  kernel  is  b\n-,k],  i.e., 

c[n;  A:]  =  ^  b[n;  Z]  a[n  —  l;k  —  Z].  (11.45) 

l 


40 


Finally,  suppose  the  c[n;  k]  and  6[n;  k]  are  both  admissibly  ergodic  kernels  in  the  sense  of 
Definition  2.  Then,  as  Vh  co  (infinite  dispersion)  and  x/i  oo  (perfect  partitioning^^ ), 
we  have,  for  each  n, 

Xm[n]^  E  [C]  Xm [n]  +  Uw  W ,  (11.46) 

where  the  Vm[n]  are  mutually  uncorrelated,  zero-mean,  complex-valued,  marginally  Gaus¬ 
sian  white  noise  sequences  that  are  uncorrelated  with  the  input  symbol  sequences 
Furthermore,  the  variance  of  the  noise  Vm[n]  is  independent  of  both  m  and  n,  and  is  given 
by 

var  Vm[n]  =  A/q >Vo£^  [l-5|^]  +  var  [C]S,  (11.47) 

where 

j  M 

fc=l 

is  the  average  transmitted  power.  In  both  (11-46)  and  (11-47)  we  have,  for  convenience, 
again  omitted  specification  of  the  time  sample  n  and  frequency  ui  corresponding  to  B  and  C 
due  to  stationarity. 


While  we  postpone  more  general  remarks  imtil  Section  II.4.2,  we  emphasize  that  The¬ 
orem  2  implies  that  the  composite  forward  link  channels  are  asymptotically  mutually  un¬ 
correlated,  identical,  nonfading,  quasi-Gaussian  chaimels  with  no  intersymbol  interference. 
Moreover,  from  Theorem  2  we  also  note  that  the  SNR  on  the  mth  of  these  composite 
channels  is 


Sm\E[C]\^ 

MomE[\B\‘^]  +  £vax[C]' 


(11.49) 


The  design  of  a  suitable  equalizer  for  the  receiver  in  this  case  is  straightforward.  In 
particular,  for  sufficiently  slow  fading  we  have 


C(w;  n]  w  A{uj-,  n]  B{uj',  n]. 


so  that  (11.49)  becomes 


^  [  \E[AB]\^ 

£  _{.hfoWo/£)E  [1B|2]  -h  var  [AB] 


(11.50) 


Then,  recognizing  that  the  term  in  brackets  in  (11.50)  is  identical  in  form  to  an  SNR  ex¬ 
pression  that  was  maximized  in  Part  I,  we  immediately  deduce  that  (11.49)  is  maximized 
when 


B{uj;n]  oa 


A*  {u;  n] 
1-l-d  ’ 


(11.51) 


where 


n] 


£\A{ui-,ri\(^ 

A7o>Vo 


(11.52) 


Actually,  perfect  partitioning  is  strictly  spealcing  not  required  for  the  forward  link  theorem  to  hold, 
although  this  is  not  apparent  in  our  proof.  However,  it  is  necessary  for  the  reverse  link  theorem  to  hold. 
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Furthermore,  in  this  case  the  corresponding  SNR  can  be  expressed  as 


(11.53) 


which  may  be  further  simplified  using  methods  used  in  Part  I.  In  fact,  as  we  will  discuss 
further  in  Section  11.5.1,  when  £m  =  £,  these  results  specialize  to  precisely  those  obtained 
in  Part  I. 


The  Reverse  Link 

The  following  theorem  summarizes  the  characteristics  of  the  composite  channel  on  the 
reverse  link.  We  emphasize  that  in  this  scenario  we  are  assuming  completely  uncoordinated 
(i.e.,  fully  asynchronous)  transmissions  from  mobiles  to  the  base.  A  proof  is  again  provided 
in  Appendix  III.D. 

Theorem  3  (Reverse  Link)  Let  XTn[n\  be  mutually  uncorrelated  sequences  of  zero-mean, 
uncorrelated  symbols,  each  with  energy  £m,  a,nd  let  am\n',k]  and  w[ri\  be  as  defined  in  (11.2) 
with  the  channel  kernels  being  statistically  independent.  Furthermore,  let  Cmiln;  fc]  denote 
the  kernel  of  the  linear  system  formed  by  cascading  the  channel  seen  by  user  m  whose 
kernel  is  amin^k],  with  the  equalizer  corresponding  to  the  detection  of  user  i,  whose  kernel 
is  hi[n-,  A:],  i.e., 

Cmi[n-,  A:]  =  ^  bi[n-,  /]  am[n  -l;k-  1].  (11.54) 

l 

Finally,  suppose  the  Cmi[n',k]  and  bmln;k]  are  admissibly  ergodic  kernels  in  the  sense  of 
Definition  2.  Then,  as  T>h  oo  (infinite  dispersion)  and  Xh  oo  (perfect  partitioning), 
we  have,  for  each  n, 

£m[n]  ^  E  [Cmm]  Xm[n]  +  «mN,  (11.55) 

where  the  Vm[^  ^*’’6  mutually  uncorrelated,  zero-mean,  complex-valued,  marginally  Gaus¬ 
sian  white  noise  sequences  that  are  uncorrelated  with  the  input  symbol  sequences  Xm\n\. 
Furthermore,  the  variance  of  the  noise  samples  Vm[ri\  is  independent  ofn  and  satisfies 

,  M 

v^Vm[n]=M^WoE[\Bm?]+  —  Y,£k^[Ckm]-  (11.56) 

Again,  in  both  (11.55)  and  (11.56)  we  have  omitted  specification  of  the  time  sample  n  and 
frequency  uj  corresponding  to  and  Ckm  due  to  stationarity. 

Again,  although  we  defer  further  more  general  remarks  until  Section  II.4.2,  we  empha¬ 
size  that,  as  in  the  forward  link  case.  Theorem  3  implies  that  the  composite  reverse  link 
channels  are  asymptotically  mutually  imcorrelated,  nonfading  quasi- Gaussian  channels  with 
no  intersymbol  interference. 
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Prom  Theorem  3,  we  observe  that  the  SNR  of  the  composite  channel  seen  by  the  mth 
user  is 

_  _ £m\E[Cmm]  I _  /-rj 

7m  —  JL4  ■  (I1.57J 


,  M 

AfoWoE  [\Bm\‘^]  +  ^  E 
k=l 


As  in  the  forward  link,  it  is  possible  to  develop  equalizers  yielding  the  best  possible  SNR 
in  the  composite  channels  for  the  reverse  link.  In  particular,  given  sufficiently  slow  fading 
that 

Cki{u}-,  n]  «  Akiuj]  n]  n], 
we  can  re-express  (11.57)  in  the  form 


Sm\E[AmBm]f 


,  M 

MoWoE  [\Bm\‘^]  +  —  '£Skvax{Ak  Bm) 

fc=i 


(11.58) 


In  turn,  the  following  lemma  directly  establishes  the  optimum  equalizer  for  the  mth  user, 
viz., 


B{u>-,n]  oc 


^  M  ’ 

1  + 


(11.59) 


fe=i 


where 


0!m(w;n]  = 


(11.60) 


A/oWo 

is  the  SNR  at  the  receiver  corresponding  to  the  mth  user.  A  proof  of  the  lemma  is  provided 
in  Appendix  III.E. 

Lemma  1  Let  Ai,A2,...,Am  denote  a  collection  of  statistically-independent,  complex¬ 
valued,  zero-mean  Gaussian  random  variables,  each  with  finite  variance,  and  let  AoWo  and 
€1,62,  ■  ■  ■  ,Sm  be  arbitrary  real,  non-negative  weights.  Then  the  function  defined 

in  (11.58)  satisfies 


1m{Brn) 

M 


< 


/ 

1  E 

\ 

E 

\ 

K 

/ 

-1 


(11.61) 


where  the  otk  are  as  defined  in  (11.60).  Furthermore,  equality  in  (11.61)  holds  when 

A* 


B  oc 


,  M 


ft=l 
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General  Remarks 

Several  aspects  of  both  the  interpretation  and  the  implications  of  Theorems  2  and  3  are 
worth  developing  in  more  detail. 

First,  it  is  important  to  emphasize  that  these  theorems  do  not  establish  that  the  compos¬ 
ite  channels  are  truly  Gaussian  channels.  For  example,  although  we  argue  that  the  additive 
noise  sequence  in  the  composite  channel  models  is  marginally  Gaussian,  we  do  not  assert 
that  these  noise  samples  are  actually  jointly  Gaussian.  More  generally,  while  these  theo¬ 
rems  describe  the  second-order  properties  of  the  composite  channels,  it  should  be  stressed 
that  they  say  very  little  about  higher-order  statistical  dependencies.^^  Indeed,  while  they 
establish  that  the  collection  of  noises  and  coded  streams  are  mutually  uncorrelated,  there 
is  no  claim  of  full  mutual  independence.  Nevertheless,  as  we  will  see  in  Section  II.5,  use¬ 
ful  approximations  are  obtained  by  modeling  the  composite  channel  as  truly  Gaussian  and 
thereby  ignoring  these  residual  statistical  dependencies. 

It  is  also  useful  to  note  that  the  additive  noise  in  the  composite  channel  model  can  be 
viewed  as  consisting  of  three  components.  Indeed,  from  both  (11.47)  and  (11.56)  we  see  that 
one  component  is  due  to  the  receiver  noise  in  the  original  fading  channel,  while  another  is 
due  to  co-channel  interference  from  other  users.  However,  the  third  component  is  a  form 
of  seif-interference  generated  as  a  by-product  of  the  modulation,  and  thus  the  total  noise 
power  has  a  term  that  depends  on  the  transmitted  signal  power.  Collectively,  these  factors, 
in  turn,  give  rise  to  the  rather  special  equalizer  structure  which  is  optimum  for  such  systems. 
In  particular,  we  note  that  the  optimum  equalizers  for  these  systems  are  minimum  mean- 
square  error  type  equalizers  in  contrast  to  the  matched  filter  t3qje  equalizers  associated  with 
traditional  CDMA  systems. 

It  should  also  be  emphasized  that  in  establishing  Theorems  2  and  3,  both  the  dispersion 
and  partitioning  characteristics  of  the  signature  set  are  important,  as  the  proofe  of  these 
theorems  reveal.  In  particular,  it  is  perfect  dispersion  that  leads  to  the  absence  of  fading  and 
intersymbol  interference  in  the  composite  channels,  while  it  is  perfect  partitioning  that  leads 
to  mutually  uncorrelated,  identical  white  noise  characteristics  in  the  composite  channels. 

In  addition,  while  Theorems  2  and  3  establish  asymptotic  results  which  involve  signature 
sequences  having  infinite  temporal  spread  (N  — >•  oo),  it  is  important  to  note  that  the 
asymptotic  behavior  can,  in  fact,  be  approximated  arbitrarily  closely  with  reahzable,  finite- 
spread  signature  sets.  In  fact,  it  this  approximation  property  that  imderlies  the  practical 
importance  of  these  theorems. 

In  general,  the  temporal  spread  required  to  achieve  a  given  level  of  approximation  to 
the  results  of  the  theorems  depends  on  the  coherence  time  characteristics  of  the  fading. 
Specifically,  analogous  to  the  results  for  spread-response  precoding,  if  N'  is  the  spread 
required  in  the  case  of  memoryless  fading,  then  the  spread  required  when  the  coherence 
time  (in  samples)  is  Ta  is 

N  =  N'(Ta  +  1). 

Thiis,  a  larger  coherence  time  implies  that  greater  delay  must  be  incurred  to  achieve  a 
given  level  of  performance.  This,  of  course,  is  true  of  communication  over  fading  channels 


fact,  because  of  the  emphasis  on  second-order  properties,  versions  of  these  theorems  hold  even  when 
the  original  channel  noises  and  kernels  are  not  Gaussian. 
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in  general.  However,  it  is  important  to  recognize  that  the  delay  constraints  in  spread- 
signature  CDMA  are  no  worse  than  in  other  systems  employing,  for  example,  interleaving. 

Although  greater  coherence  times  require  longer  signature  sequences,  it  is  important  to 
note  that  the  computational  requirements  need  not  grow  with  coherence  time.  In  fact,  the 
required  number  of  non-zero  coefficients  in  each  signature  sequence  is  independent  of  the 
coherence  time.  In  particular,  to  match  the  coherence  time  characteristics  of  a  channel  of 
interest,  it  suffices  to  upsample  some  prototype  set  of  signature  sequences  (such  as  those 
developed  in  Section  II.3.1)  to  achieve  the  necessary  temporal  spread.  Furthermore,  we 
note  that  the  energy,  dispersion,  partitioning,  and  computational  characteristics  are  all 
unaffected  by  such  upsampling. 

While  this  upsampling  approach  to  coherence  time  matching  is  also  used  with  spread- 
response  precoding,  an  important  difference  in  the  case  of  spread-signature  CDMA  is  that 
the  upsampling  factor  K  cannot  be  freely  chosen.  In  particular,  in  order  to  preserve  or¬ 
thogonality  among  the  signatures,  we  may  only  upsample  the  signatmes  by  an  integer  K 
that  is  not  a  multiple  of  a  prime  factor  of  the  number  of  users  Hence,  if  the  hm[n]  are 
a  set  of  prototype  maximally-spread  signatme  sequences  applicable  to  memoryless  fading, 
the  signature  sequences  applicable  when  the  coherence  time  is  Tq  are  given  by 


I  hm[n/Ko]  n  =  ■  •  •  ,  -Kq,  0,  Kq,  2Kq,  ■  ■  • 
1  0  otherwise 


(11.62) 


where  Kq  denotes  the  smallest  integer  that  is  not  less  than  Tq  +  1  and  is  not  a  multiple  of 
a  prime  factor  of  M. 


II.  5  Performance 

In  this  section  we  explore  both  the  effective  capacity  and  bit-error  rate  characteristics 
of  spread-signature  CDMA  systems.  Again,  we  examine  the  forward  and  reverse  links 
separately  so  as  to  emphasize  some  important  differences  between  these  cases. 

II.5.1  The  Forward  Link 

Due  to  the  coordinated  nature  of  the  transmission,  performance  on  the  forward  link  gener¬ 
ally  follows  immediately  from  the  results  presented  for  single-user  systems  in  Part  I.  Indeed, 
on  the  forward  link,  it  is  typical  to  choose 

Si=€2  =  ---  =  £m  =  €  =  £, 

in  which  case  the  optimum  equalizer  (11.51)  and  the  optimum  SNR  (11.53)  both  specialize 
to  precisely  those  derived  in  the  single-user  scenario.  Furthermore,  since  forward  link  per¬ 
formance  for  any  M  coincides  with  reverse  link  performance  for  M  =  1,  we  will  ultimately 
rederive  these  results  as  special  cases  in  Section  II.5.2. 


^®This  implies,  for  instance,  that  when  the  number  of  users  M  is  a  power  of  two,  i.e.,  M  =  2"^,  we  can  use 
any  odd  K.  Likewise,  for  M  =  ,  we  can  use  any  K  which  is  not  a  multiple  of  3,  and  for  M  =  ,  we  can 

use  any  K  which  is  not  a  multiple  of  2  or  3,  i.e.,  K  =  5,7, 11, . . .,  etc. 
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Figure  11.4:  Bit  error  probabilities  using  uncoded  QPSK  in  Rayleigh  fading  on  the  for¬ 
ward  link.  The  dash-dotted  and  dashed  curves  represent  the  performance  obtained  using 
maximally-spread  signature  sets  with  N'  =  32  and  N'  =  128,  respectively.  The  solid  line 
represents  the  ideal  performance  bound  (N'  —>  oo). 


While  many  of  the  performance  characteristics  on  the  forward  link  may  be  directly  ex¬ 
tracted  from  the  corresponding  results  in  Part  I,  those  concerning  performance  with  realiz¬ 
able  (finite-spread)  signature  sets  cannot.  Accordingly,  we  depict  forward  Unk  performance 
obtained  with  realizable  signature  sets  in  Fig.  II.4.  Specifically,  we  plot  realizable  uncoded 
quadrature  phase-shift  keying  (QPSK)  bit  error  rate  performance  in  Rayleigh  fading  as  a 
function  of  SNR  per  bit  as  measured  through  Monte  Carlo  simulations.  In  these  simulations 
we  use  the  maximally-spread  signature  sets  of  dimension  M  —  2  developed  in  Section  II.3.1, 
and  the  dash-dotted  and  dashed  curves  indicate  the  performance  corresponding  to  N'  =  32 
and  N'  =  128,  respectively.  For  comparison,  the  solid  curve  corresponds  to  the  asymptotic 
performance  obtained  with  infinite  temporal  spread  {N'  oo). 

More  generally,  the  single-user  scenario  (M  =  1)  warrants  further  discussion.  As  re¬ 
marked  earlier,  no  nontrivial  finite-spread  orthogonal  signatures  for  M  =  1  exist.  As  a 
result,  in  Part  I  the  orthogonality  (losslessness)  constraint  had  to  be  relaxed  in  develop¬ 
ing  realizable  precoders.  In  turn,  this  led  to  a  corresponding  degradation  in  the  precoder 
performance.  However,  from  within  the  broader  framework  of  this  report,  we  can  see  that 
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orthogonal  signature  sets  corresponding  to  M  >  2  can  also  be  usefully  applied  to  single-user 
precoding  problems.  In  particular,  we  may  distribute  successive  symbols  in  a  stream  x[n] 
among  M  different  substreams  Xm[n]  by  periodic  subsampling,  ?.e., 

a^mN  =  x\nM  ■+■  m  —  1]. 

and  then  xise  an  orthogonal  spread  signature  set  of  order  M  to  implement  perfectly  loss¬ 
less  precoding.  At  the  receiver,  the  substreams  may  be  demodulated  from  the  signature 
sequences,  and  appropriately  reintegrated  into  the  single  stream  x[n]. 

This  strategy  for  implementing  lossless  precoding  has  a  conceptually  rather  useful  in¬ 
terpretation.  Specifically,  the  cyclic  distribution  of  symbols  among  substreams  may,  itself, 
be  interpreted  as  a  lossless  system,  and,  hence,  in  effect  what  we  construct  in  this  manner 
is  a  linear  periodically  time- varying  precoder Hence,  we  see  that  by  adopting  this  more 
general  precoder  structure  we  are  able  to  achieve  lossless  spread-response  precoding  with 
finite-length  filters. 

In  light  of  these  observations,  we  see  that  the  curves  in  Fig.  II.4  also  correspond  to  the 
single-user  performance  that  is  achievable  using  the  improved  precoders.  And,  comparing 
Fig.  II.4  to  Fig.  1.5,  we  see  that  significantly  better  performance  is  achieved  for  a  given 
length  constraint  (iV')  when  we  have  perfect  orthogonality  (losslessness). 

II.5.2  The  Reverse  Link 

The  performance  characteristics  on  the  reverse  link  are  dramatically  different  than  those  on 
the  forward  link  due  the  the  uncoordinated  nature  of  the  transmissions,  as  we  now  illustrate. 
We  begin  by  noting  that,  in  practice,  on  the  reverse  link  power  control  is  generally  employed 
to  adjust  the  transmitter  power  for  each  user  so  that  the  average  received  power  from  each 
user  is  the  same,  i.e.,  for  all  users  m,  the  SNR  in  the  corresponding  original  channel  is 

E[am{oj;n]]  =  l/(:o.  (11.63) 

where  am(<^;7^]  is  as  defined  in  (11.60).  Indeed,  such  power  control  is  a  practical  technique 
for  mitigating  the  vulnerability  of  CDMA  systems  to,  among  other  problems,  near-far  effects 
[18]. 

Some  important  insights  into  the  potential  capacity  of  spread-signature  CDMA  systems 
arise  out  of  the  theorems  of  Section  II.4.2.  Theorem  3,  in  particular,  suggests  that  the 
composite  quasi-Gaussian  channel  seen  by  the  mth  user  has  a  “capacity”  given  by 

Cm  =  Wo  log(l  4-  7m)  (11.64) 

where  •jm  is  the  SNR  in  the  composite  channel  of  the  mth  user.  However,  we  emphasize 
that  caution  must  be  exercised  in  interpreting  this  notion  of  capacity.  In  particular,  we  note 
that  this  measure  of  capacity  ignores  aU  statistical  effects  of  greater  than  second  order  in 
the  composite  channels — despite  the  fact  that  such  effects  can  in  principle  be  exploited  to 


^■^We  remark  as  an  aside  that  traditional  interleavers  are  essentially  trivial  examples  of  such  linear  peri¬ 
odically  time-varying  systems. 
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yield  still  higher  throughputs.  Nevertheless,  the  particular  concept  of  constrained  capacity 
defined  here  has  a  useful  interpretation  as  we  will  develop.  In  particular,  it  provides  a  mea¬ 
sure  of  throughput  that  can  be  achieved  when  traditional  forms  of  coding  are  superimposed 
on  our  CDMA  system. 

Convenient  expressions  can  be  obtained  for  the  composite  channel  SNR  which  maximizes 
this  effective  capacity.  In  particular,  from  Lemma  1,  we  have  that  when  we  use  optimum 
equalization,  7^  can  be  expressed  as 


7m  _  I 

'M~  Pm 


(11.65) 


where 


(11.66) 


Under  our  power  control  assumption.  Pm  as  defined  in  (11.66)  is  independent  of  m,  and  thus 
we  simply  use  Pq,  70,  and  C  to  denote  the  quantities  (11.65),  (11.66),  and  (11.64)  respectively. 
Furthermore,  in  this  case,  by  exploiting  symmetry  (11.66)  can  be  expressed  in  closed  form 
in  terms  of  the  exponential  integral.  In  particular,  we  have  the  following  Lemma,  whose 
proof  is  provided  in  Appendix  III.F. 


Lemma  2  Let  vi,V2,  ■  ■ .  ,vm  be  independent,  identically-distributed  exponential  random 
variables  with  mean  l/p.  Then  for  any  1  <m  <  M, 


k^m  _  M-1 
1  +  M  Ml 

k 


r  ^“2  . ,  ] 

L  fc=o  ^  J 


(11.67) 


where 

fOO  g-t 

JSi(i/)=/  —df.  (11.68) 

Ji/  t 

is  the  exponential  integral  [15]. 

Thus,  using  Lemma  2  with  Vk  =  we  immediately  obtain  that  (11.66)  can  be 

expressed  as 


,  M-l  ,  (VCo)" 


(-l)^+^e^^'>Ui(MCo)  +  Y.  (-1)^ 


k\ 

(MCo)^+i 


(11.69) 


Hence,  evaluating  (11.64)  via  (11.65)  and  (11.69),  we  obtain  relationships  for  the  capacity 
per  user  as  a  function  of  the  number  of  users  M,  the  SNR,  and  the  bandwidth  per  user. 
In  Fig.  11.5  we  plot  capacity  per  user  as  a  function  of  M  at  various  levels  of  SNR,  while  in 
Fig.  II.6  we  plot  capacity  per  user  as  a  function  of  SNR  for  various  values  of  M,  with  the 
Gaussian  channel  capacity  included  for  comparison. 

Note  firom  Fig.  II.6  that  for  M  >  2  the  reverse  link  capacity  appears  to  saturate  at 
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Number  of  Users 


Figure  II.5;  Reverse  link  capacity  per  user  as  a  function  of  the  number  of  users  M.  The 
successively  higher  curves  correspond  to  SNRs  of  -5,  5,  15,  and  25  dB.  Note  that  the 
capacities  corr^ponding  to  M  =  1  also  coincide  with  forward  link  performance  with  any 
number  of  users.  The  connecting  dashed  lines  have  no  special  interpretation;  they  are 
provided  as  visual  aides  only. 
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Figure  II.6:  Capacity  per  user  as  a  function  of  SNR.  The  successively  lower  solid  curves 
correspond  to  reverse  link  transmission  with  M  =  1,  M  =  2,  M  =  A,  and  M  — )•  oo  users. 
Note  that  the  curve  for  M  =  1  also  coincides  with  forward  link  capacity  with  any  number 
of  users.  The  dashed  line  indicates  the  capacity  of  the  corresponding  Gaussian  channel. 
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high  SNR,  with  the  saturation  level  depending  on  the  number  of  users.  This  is,  in  fact, 
the  case,  and  it  reflects  the  fact  that  reverse  link  performance  is  fundamentally  interference 
limited  rather  than  noise  limited.  This  phenomenon,  which  is  a  direct  consequence  of  the 
constrained  receiver  structure  we  employ,  is  well  known;  see,  e.g.,  Verdu  [33].  While  it 
is  possible  to  design  multiple-access  systems  that  are  strictly  noise  Mmited,  this  requires 
the  use  of  receivers  employing  joint  detection  strategies.  Unfortunately,  joint  detection  is, 
in  general,  computationally  very  intensive,  particularly  when  the  number  of  users  in  the 
system  is  large  [34].  For  these  reasons,  such  receivers  are  often  considered  impractical. 

To  verify  the  interference-limited  behavior,  let  us  consider  the  high  SNR  regime  perfor¬ 
mance  in  more  detail.  In  this  case,  we  may  use  the  series  expansion  [15] 


OO 


Ei(i/)  = -Fo  -  Ini^  - 

fc=i 


(-1)^1/* 

kk\ 


(11.70) 


with  To  =  0.57721  •  •  •  denoting  Euler’s  constant,  to  show  that  at  high  SNR  (Co  — t  0),  (11.69) 
satisfies 

M-  1 

/3o  (11.71) 

Thus  the  asymptotic  capacity  on  the  reverse  link  in  this  regime  is 

Hence,  (11.72)  implies  that  the  resulting  capacity  is  finite,  consistent  with  the  saturation 
behavior  apparent  in  Fig.  11.6. 

Another  important  asymptotic  regime  to  explore  is  the  large  number  of  user  scenario. 
Fortunately,  as  M  -)•  oo,  convenient  closed  form  expressions  are  possible.  In  particular, 
since  the  ak  are  independent,  identically-distributed  random  variables,  we  have,  by  the  law 
of  large  numbers 

1  ^ 

VCo  as  M  OO. 

k=\ 

Hence,  we  get 

7.  ^  (U.73) 

and,  in  tirm,  (11.64)  becomes 

C -4  Wo  log  (11.74) 

Note  from  (11.74)  that  when  the  system  contains  large  nmnbers  of  users,  we  have,  regardless 
of  the  transmitter  power  used,  C/Wo  <  1  bit/sec/Hz  with  C/Wo  1  bit/sec/Hz  as  Co  0. 

It  is  also  significant  that  in  this  case,  the  equalizer  structure  simplifies  substantially  as 
well.  Again,  via  the  law  of  large  numbers  we  get 


R  oc  A 


♦ 


(11.75) 


the  famihar  matched  filter  equalizer,  which  we  note  is  substantially  easier  to  implement  than 


51 


Figiire  II.7:  Capacity  per  user  as  a  function  of  bandwidth.  The  successively  lower  solid 
curves  correspond  to  reverse  link  transmission  with  M  =  1,  M  =  2,  M  =  4,  and  M  —>■  oo 
users,  respectively.  Note  that  the  M  =  1  curve  for  the  reverse  link  coincides  the  performance 
of  the  forward  link  with  any  number  of  users.  The  dashed  line  corresponds  to  the  capacity 
of  the  Gaussian  channel. 

the  optimum  equalizer  for  finite  M.  Fortunately,  the  asymptotic  performance  is  achieved 
for  fairly  moderate  values  of  M  as  Fig.  11.5  reveals,  so  that  in  practice  this  simpler  equalizer 
can  frequently  be  used. 

Let  us  now  turn  our  attention  to  the  SNR-limited  setting,  and  consider  the  large  band¬ 
width  behavior.  In  Fig.  II.  7  we  plot  normalized  capacity  as  a  function  of  bandwidth  per 
user  for  various  values  of  M,  with  the  Gaussian  channel  capacity  included  for  comparison. 
Several  general  remarks  regarding  Fig.  II.7  are  appropriate.  First,  as  in  the  case  of  the 
Gaussian  channel,  the  capacity  in  our  multiuser  fading  environment  increases  monotoni- 
cally  with  bandwidth  per  user.  Second,  the  horizontal  distance  between  the  M  =  1  and 
M  >2  curves  quantifies  how  much  additional  bandwidth  is  required  to  compensate  for  the 
absence  of  coordination  on  the  reverse  link. 

Perhaps  the  most  important  observation  from  Fig.  II.  7,  however,  is  that  the  infinite 
bandwidth  capacity  is  identical  to  that  of  the  Gaussian  channel.  To  verify  this  property, 
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we  note  that  using  the  asymptotic  series  expansion  [15] 


k\ 


Ei{v)  =  e  large  u, 


(11.76) 


fe=0 


we  have  that  for  large  Wq  {i.e.,  large  Co)? 


fc=o  (^Co)''+i' 


(11.77) 


Using  (11.77)  in  (11.65),  we  see  that  for  large  Wq, 

M 


70 


MCo-1 

However,  since  70  in  (11.78)  is  small,  (11.64)  can  be  expressed  asymptotically  as 

MWo 


C  ~  Wo7o  = 


Using  the  fact  that  (11.60)  and  (11.63)  imply 


MCo-1' 


(11.78) 


(11.79) 


we  verify  from  (11.79)  that 


IV  jc  —  [i^lmP] 
“ - Me - 


^^0 


(11.80) 


which  we  emphasize  is  independent  of  m  by  our  power  control  assumption.  Finally,  we  note 
that  (11.80)  is,  of  course,  the  capacity  of  the  infinite-bandwidth  Gaussian  channel. 


II.5.3  Exploiting  Additional  Processing  Gain 

For  reverse  link  communication,  co-channel  interference  is  such  that  even  at  high  SNR,  the 
SNR  7o  of  the  composite  quasi-Gaussian  channel  for  each  user  is  low.  In  particular,  using 
(11.71)  in  (11.65)  we  see  that  regardless  of  the  available  bandwidth  or  power  we  have 


0  <  7o  < 


M 

M-1 


for  M  >  2.  This  means  that  in  order  to  achieve  a  throughput  per  user  anywhere  near  that 
predicted  by  the  capacity  calculations  in  Section  II.5.2,  an  enormous  amount  of  coding  is 
required.  Conversely,  without  the  use  of  coding,  bit  error  rate  performance  on  the  reverse 
link  will  be  invariably  poor. 

As  an  alternative  to  coding,  one  can  consider  exploiting  simple  spread-spectrum  process¬ 
ing  gain  to  raise  70  to  a  level  sufficient  for  acceptable  bit  error  rate  performance.  Although, 
as  we  shall  see,  direct  coding  is  always  preferable  in  terms  of  efficiency,  using  simple  pro¬ 
cessing  gain  is  appealing  in  terms  of  its  substantial  ease  of  implementation  and  reduced 


53 


computational  complexity. 

Augmenting  our  M-user  system  model  to  incorporate  such  processing  gain  is  straight¬ 
forward.  In  essence,  we  expand  the  bandwidth  per  user  while  fixing  the  symbol  rate  by 
replacing  the  up-  and  down-sampling  rates  in  Figs.  II.  1  and  II.3  with  an  integer  L  satisfying 
L  >  M.  The  resulting  bandwidth  expansion  factor  or  processing  gain  is  then 

p  =  L/M.  (11.81) 


With  this  change,  the  SNR  7m  of  the  composite  channel  becomes 

7m  _  I  _  , 

L~  13m  ’ 


where 


Pm  =  E 


■i  m 


,  M 

l  +  jf^ak 


(11.82) 


and  with  ak  as  in  (11.60).  Again,  using  Lemma  2  with  Vk  =  cx-k/L  we  can  express  (11.82)  in 
closed  form  as 

In  the  remainder  of  this  section,  we  focus  on  the  large  nmnber  of  user  (M  ->•  oo)  scenario. 
In  this  case,  again  using  the  law  of  large  numbers,  we  get 


Co  +  l/p’ 


(11.84) 


where  p  is  as  given  by  (11.81).  At  high  SNR  {i.e.,  when  the  receiver  noise  is  negligible), 
(11.84)  further  simplifies  to 

7o  ~  P- 

The  use  of  processing  gain  permits  a  tradeoff  between  bit  error  rate  and  bandwidth.  In 
particular,  from  (11.84)  we  see  that  larger  processing  gains  invariably  give  rise  to  improved 
bit  error  rate  performance.  However,  this  improvement  is  obtained  at  the  expense  of  band¬ 
width.  Prom  an  information  theoretic  point  of  view,  this  particular  tradeoff  is  not  efficient. 
Too  see  this,  note  that  if  we  apply  coding  on  top  of  this  new  system  that  exploits  processing 
gain,  the  apparent  capacity  per  user  per  unit  bandwidth  is 


=  llog(l+7.)  =  bog(l  +  ^^) 


(11.85) 


where  the  last  equality  follows  from  (11.84).  In  Fig.  II.  8,  we  use  (11.85)  to  plot  capacity  per 
user  per  unit  bandwidth  as  a  function  of  the  processing  gain  p  in  the  large  number  of  users 
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Figure  11.8:  Capacity  per  user  per  unit  bandwidth  as  a  function  of  processing  gain  p,  where 
the  number  of  users  satisfies  M  — >  oo.  The  successively  higher  curves  correspond  to  SNRs 
of  -5,  5,  15,  and  25  dB. 

limit  M  oo.  We  note  that  the  effective  capacity  achievable  with  additional  coding  falls 
monotonically  with  p,  and  that  the  case  p  =  1  corresponds  to  our  original  system,  i.e.,  the 
system  without  processing  gain.  This  is  because  the  bandwidth  expansion  incurred  by  using 
a  larger  processing  gain  p  more  than  offsets  the  increase  in  the  SNR  per  symbol  achieved. 
In  fact,  it  is  straightforward  to  show  from  (11.85)  that  C  — >  0  as  p  oo. 

Despite  the  inherent  bandwidth  inefficiency,  exploiting  processing  gain  is  highly  prac¬ 
tical  technique  and  enjoys  widespread  use.  Let  us  consider,  for  example,  the  bit  error 
rate  performance  of  uncoded  QPSK  (quadrature  phase-shift  keying)  using  spread-signature 
CDMA  with  M  users  and  processing  gain  p  =  L/M,  some  integer  L.  In  this  case,  again 
assuming  the  use  of  power  control,  the  bit  error  probability  is  given  by 

V=Q{VjE)  (11.86) 

where 
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and  where  70  is  given  by  (II.5.3)  with  (11.83),  or,  in  the  case  of  large  numbers  of  users, 
(11.84). 

In  Fig.  II.  9  we  plot  bit  error  probability  as  a  function  SNR  per  bit  for  several  processing 
gains  in  the  large  number  of  user  scenario,  and  include  comparisons  to  conventional  CDMA. 
In  both  the  spread-signature  and  conventional  CDMA  systems,  no  channel  coding  is  used. 
Prom  this  plot,  it  is  apparent  that  there  is  an  enormous  advantage  in  using  spread-signature 
CDMA  over  conventional  CDMA  in  such  uncoded  systems.  This  is  due  to  the  fact  that 
spread-signature  CDMA  is  much  more  effective  than  conventional  CDMA  at  mitigating  the 
effects  of  fading  even  without  additional  coding.  In  particular,  at  high  SNR  and  for  large 
numbers  of  users,  the  bit  error  probability  of  spread-signatme  CDMA  saturates  at 

^~Q(Vp),  (11-87) 

while  for  conventional  CDMA  the  bit  error  probability  saturates  at 


Finally,  in  Fig.  II.  10  we  plot  bit  error  probability  as  a  function  of  processing  gain  for 
several  SNR  values,  again  in  the  large  number  of  users  scenario.  For  comparison,  we  also 
plot  the  performance  of  conventional  CDMA,  using  (11.88). 

II.6  Summary 

Spread-signature  CDMA  as  developed  in  this  part  of  the  report  constitutes  a  potentially 
attractive  alternative  to  conventional  CDMA  for  general  multiuser  communication  in  fad¬ 
ing  environments.  Indeed,  the  performance  results  suggest  that  by  effectively  combining 
modulation  and  precoding,  such  systems  appear  to  offer  significantly  better  performance. 

Nevertheless  numerous  issues  remain  to  be  explored.  In  traditional  CDMA  systems, 
coding  is  used  to  combat  the  effects  of  fading,  co-channel  interference,  and  receiver  noise. 
However,  in  spread-signature  CDMA  systems,  the  effects  of  fading  are  mitigated  by  the 
precoding  implicit  in  the  modulation,  leaving  coding  to  handle  only  the  remaining  interfer¬ 
ence  and  noise  effects.  This  would  appear  to  be  a  favorable  computational  tradeoff,  since 
in  general  demodulation  is  computationally  much  cheaper  than  decoding.  However,  further 
simulations  are  required  to  verify  these  computational  advantages.  In  addition,  while  we 
have  restricted  our  attention  to  single-user  detector  receivers,  the  possibility  of  developing 
viable  joint  detector  receivers  is  also  worth  pursuing. 

Other  issues,  both  technical  and  practical,  also  warrant  further  investigation.  Among 
the  technical  issues,  an  obvious  question  that  arises  concerns  the  extent  to  which  the  capac¬ 
ity  expressions  developed  in  this  report  can  be  interpreted  in  a  strict  Shannon  sense.  Some 
of  the  more  important  practical  issues  which  remain  to  be  explored  include  receiver  S3m- 
chronization  and  peak-to-average  transmitter  power  requirements,  and  equalizer  sensitivity 
characteristics. 
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Figure  11.9:  Bit  error  probability  as  a  function  of  SNR  per  bit  for  uncoded  QPSK  on  reverse 
link  with  M  oo  users.  The  successively  lower  solid  curves  correspond  to  the  performance 
of  spread-signature  CDMA  with  processing  gains  of  p  —  1, 7, 13, 19.  For  comparison,  the 
successively  lower  dashed  curves  correspond  to  the  performance  of  conventional  CDMA 
with  the  same  series  of  processing  gains.  For  these  comparisons,  no  coding  is  used  in  either 
system. 
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Figtire  II.  10;  Bit  error  probability  as  a  function  of  processing  gain  for  uncoded  QPSK  on 
reverse  link  with  M  oo  users.  The  successively  lower  solid  curves  correspond  to  the 
performance  of  spread-signature  CDMA  with  SNRs  of  -5,  5,  15,  and  25  dB/bit.  For  com¬ 
parison,  the  successively  lower  dashed  curves  correspond  to  the  performance  of  conventional 
CDMA  with  the  same  series  of  SNRs.  For  these  comparisons,  no  coding  is  used  in  either 
system. 
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Part  III 


Appendices 


III.  A  Proof  of  Theorem  1 

As  an  intermediate  step,  we  obtain  results  concerning  the  following  related  system.  Let  S  {•} 
denote  a  linear  system  which  is  the  cascade  of  an  LTI  system  whose  unit-sample  response 
pi[n],  followed  by  a  linear  system  whose  response  at  time  n  to  an  unit-sample  at  time  n  —  k 
is  f[n‘,k],  followed  by  another  LTI  system  whose  unit  sample  response  is  5'2[«]-  Hence, 

q[n]  =  S  {p[n]}  =  -  i]  -i-k-l]  (III.l) 

i  k  l 

Furthermore,  let  «[n;  fc]  denote  the  kernel  of  the  overall  linear  system,  i.e.,  its  response  at 
time  n  to  a  -unit-sample  at  time  n  —  k. 

We  begin  with  a  useful  lemma  regarding  such  systems. 


Lemma  3  Let  pi[n]  and  52  M  be  lossless^  and  let  f[n\k]  be  the  kernel  of  an  admissibly 
ergodic  system.  Then  as 


(IIL2) 

the  kernel  u[n-,  A;]  defined  above  obeys 

u[n;  k]^  fi  (pi[A:]  *  p2[fc])  =  u[A:] 

(IIL3) 

and 

u[n; n  —  k]u* [m; m  —  k]  a^5\n  —  m] 

k 

(IIL4) 

where 

ii[n;  fc]  =  u[n;  k]  —  u[A:]. 

Proof: 

Prom  (III.l)  we  get 

n[n;  k]=Y^ 52W  XI  -  i  -  1]- 

(IIL5) 

i  I 
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Then 


E  [u[n;  A:]]  =  ^ 52 W  ffi  [k-i]  =  iic-  9\[k]  * 


u[n;  fe]  =  92[i\  f[n  -i;l]gi[k-i-  1] 

i  I 


where,  consistent  with  (1.14), 


f[n;  As]  =  /[n;  k]-  fx  J[A:]. 


Hence, 


[|tt [n;  A;]  1^  =Yl'^92[i]92 [*'] [A;  - /  -  *] 51  ^  -  *']  (111.6) 


where  i2[m;A:]  is  as  defined  in  (L17b).  Applying,  in  order,  the  triangle  inequality  and  the 
Cauchy  inequality,  we  are  able  to  bound  (III.  6)  by 


where 


provided 


E  pij[n,  A:]|^]  <  <f{9i)<p{92) 


<P^i9)  =  ^  ^11  9\i']  <  Sr 

i,i* 


(III.7) 


^  1- 


Applying  the  Cauchy  inequality  again  to  <^^(52)}  however,  gives,  after  some  simplification 
and  using  (III.  2), 

^"(52)  <  Sr/V2.  (III.8) 

Hence,  (IIL8)  and,  in  turn,  (IIL6)  tend  to  zero  as  T>2  00,  which  verifies  (III.3). 

To  show  (III.4),  we  begin  by  noting  that  since  pi[n]  is  lossless,  we  get,  using  (1.15), 

p[n,  m]  =  ^  u[n\  n  —  A:]  u*[m;  m  —  A:]  =  ^  92\n  —  i]  92\fn  —  i  + 1]  d[z;  /]. 

k  i,l 


Then,  since  by  (I.17b) 


we  get,  using  (1.25), 


E  [c?[n;  A:]]  =  (5[A:], 


E  [p[n,  m]]  =  5[n  —  m]. 


Next,  we  write,  using  (1.16), 


/5[n,  m]  =  p[n,  m]  —  a‘^5[n  —  m]  =  ^  j2[«  —  *]  92['m  —  i  4-  /]  d[i;  /] 
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and,  using  (I.17c) 

E  [lp[n,  m]|^]  =  -  i]  92[m  -i  +  l]  g2[n  -  i']  g2[m  -i'  + 1']  T[i  -  i'-,  I,  l'].  (III.9) 

However,  applying  the  triangle  and  Cauchy  inequalities,  in  order,  and  noting 

T[m-,kJ]=T*[-m;l,k] 

we  can  bound  (III.9)  by 

E[\p[n,m]\^]  <  +  (IILIO) 

Applying  the  triangle  and  Cauchy  inequalities,  in  order,  again,  to  the  right  side  of  (III.  10), 
and  changing  variables,  we  get 

E  [lp[n,m]l2]  <  5IW  \T\i'-,lA\  =  St/V2.  (III.ll) 

which  approaches  zero  as  V2  00.  I 

Using  Lemma  3,  we  obtain  the  following  theorem: 

Theorem  4  Let  51  [n]  and  52  W  be  lossless,  and  let  f[n;k]  be  an  admissibly  ergodic  system 
kernel.  Further,  and  suppose  p[n]  is  a  wide-sense  stationary,  white  random  process  with 
mean  zero  and  variance  a^,  and  that  p[n]  and  the  channel  kernel  f[n;  /:]  are  statistically 
independent.  Then,  as  V2  — >  00,  g[n]  as  defined  by  (III.l)  satisfies,  with  u[n]  as  defined  in 

(III.S), 

q[n]  «[n]  *p[n]  +  z{n],  (III.12) 

where  2:[n]  is  a  marginally  Gaussian  wide-sense  stationary,  white  random  process  with  mean 
zero  and  variance  and  where  z[n]  is  uncorrelated  withp[n]. 


where  we  use  Ep  [•]  to  denote  expectation  with  respect  to  p[n]  given  a  fixed  but  arbitrary 
realization  of  the  kernel  /[n;  k]. 
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In  addition,  from  (111.13)  we  get 

Ep  [.2;[n]p*[m]]  =  ^«[n;  A:]  Ep  [p[n  —  A;jp*[m]]  =  o^u\n-,n  —  m],  (III.14) 

k 

Hence,  since  (III.  3)  in  Lemma  3  implies 

«[n;  A:]  0, 

we  get,  from  (III.  14), 

Ep  [2:[n]p*[m]]  0  as  I>2  00. 

Finally,  using  (III.  13)  we  also  obtain 

Ep  [z[n]  z* [m]]  =  ^  u[n;  A:]  u*  [m;  l]Ep  [p[n  —  A:]  p* [m  —  Z]] 
k,l 

=  a^^u[n;n  —  A:]u*[m;m  —  A:].  (III.15) 

k 

Applying  (III.4)  in  Lemma  3  to  (111.15),  we  see  immediately 

Ep  [z[n\  z*[m]\  a^o^5[n  —  m]  as  X>2  — >■  oo. 

Finally,  to  show  that  z[n]  is  marginally  Gaussian  requires  a  straightforward  Central 
Limit  Theorem  argument.  ■ 

The  following  pair  of  corollaries  lead  directly  to  a  proof  of  our  main  theorem. 


Corollary  1  Suppose  51  [n]  =  52  [—J^]  in  Theorem  4-  Then  as  V2  00, 

g[n]  Aip[n]  +  z\n\.  (111.16) 


Proof: 

It  suffices  to  note  that,  since  g2[n]  is  lossless 

51 N  *  52 W  =  52 W  *  52[-w]  = 
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Furthermore,  since  52 N  is  lossless,  p[n]  and  p[n]  *  P2W  have  the  same  spectrum  and  thus 

varp  (PP2N  *p[n])  =  \fi\‘^cr^. 

Applying  (1.25)  and  (1.19)  we  then  verify  (III.17).  ■ 

Finally,  to  establish  Theorem  1  we  need  only  recognize  that  x[n]  can  be  divided  into 
two  uncorrelated  components; 

x[n]  =  xi[n]  +  X2[n], 

where  Xi[n]  is  the  component  due  to  x[n]  and  where  X2[n]  is  due  to  u;[n].  Using  Theorem  4 
and  Corollary  1  with  p[n]  =  x[n],  9[n]  =  xi[n],  and  f[n;  A:]  =  c[n;  fc],  we  see 

xi[n]  +  z\n\ 

where  the  white  noise  z[n]  has  variance  SsO^.  Similarly,  using  Theorem  4  and  Corollary  2 
with  p[n]  =  'w[n],  9[n]  =  X2[n],  and  f[n;  fe]  =  6[n;  k]  we  get  that 

var;r2[n]  =  AfoWoE  pB(a;;n]p]  . 


III.B  Linear  Randomly  Time- Varying  Systems 

We  adopt  the  following  notation  for  linear  randomly  time-varying  systems  corresponding 
to  wide-sense  stationary  uncorrelated  scattering.  We  begin  by  using  /[n;A:]  to  denote  the 
kernel  of  a  discrete-time  linear  system,  corresponding  to  the  response  of  the  system  at  time 
n  to  a  unit  sample  at  time  n  —  k.  Hence,  the  response  of  the  system  to  an  input  x[n]  is 

y  W  =  H  /[”;  -  ^1- 

k 

The  time- variant  system  frequency  response  associated  with  this  system  is  denoted  by 

F(a;;n]  =  ^/[n;A:]e-^■‘^^  (III.19) 

k 

and  represents  the  response  of  the  system  to  complex  exponential  For  stationary 

systems,  we  define  the  system  correlation  function  by 

R[m;  k]  =  E  |^/[n;  k]  f*[n  —  m;  A;]j  . 

where 

/[n;  k]  =  f[n;  k]-E  [/[n;  A:]] , 
and  the  system  scattering  function  by 

5(A;  k]=Y^  -RN;  k] 
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When  the  system  is  chareicterized  by  uncorrelated  scattering  we  have  both 


£J[/[n;A:]]  =  E[/[n;0]]  5[k] 


and 

E  j^/[n;  fe]  f*[n  —  m;  Z]j  =  R[m\  k]  (5[fc  —  Z]. 

This  makes  Z?(a;;n]  defined  in  (111.19)  wide-sense  stationary  in  both  uj  and  n,  a  property 
we  exploit  extensively. 

Another  useful  characterization  of  such  systems  is  in  terms  of  the  associated  spaced- 
frequency  spaced-time  correlation  function,  which  is  given  by 

m]  =  R[m\  A:]  =  E  ^{8\  n]  F* {6  —  u:-^n  —  m 

k 


where 

F(u;;  n]  =  F{u-,  n]-E  [F(a;;  n]] . 

In  turn,  we  define  the  system  Doppler  power  spectrum  by 

T(A)  =  5^«r(0;m]e^^"*. 

m 

Finally,  the  multipath  intensity  profile  or  delay  power  spectrum  of  the  system  is 

n[A:]  =  i?[0;  A;]  =  var  /[n;  A:], 


so  the  total  power  is 

varF(a;;n]  =  ^(0;0]  =  5^n[A;]. 

k 

III.C  Maximally  Spread  Signature  Sequences 

Coefficients  are  first  converted  to  a  binary  sequence,  with  1/\/N  represented  by  T’  and 
—\/\/N  represented  by  ’O’,  then  replaced  with  the  hexidecimal  equivalent. 
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M  N  h[n] _ 

2  4  E 

D 

4  4  F 

A 
C 
9 

2  8  ED 

E2 

8  8  FF 

AA 
CC 
99 
FO 
A5 
C3 
96 

2  16  EDE2 

EDID 

4  16  FAC9 

F5C6 
FA36 
F539 

2  32  EDE2ED1D 

EDE212E2 

2  64  EDE2ED1DEDE212E2 

EDE2ED1D121DED1D 

4  64  FAC9F5C6FA36F539 

FAC90A39FA360AC6 
FAC9F5C605C90AC6 
FAC90A3905C9F539 

8  64  FFAACC99F0A5C396 

FF55CC66F05AC369 
FFAA3366F0A53C69 
FF553399F05A3C96 
FFAACC990F5A3C69 
FF55CC660FA53C96 
FFAA33660F5AC396 
FF5533990FA5C369 

2  128  EDE2ED1DEDE212E2EDE2ED1D121DED1D 

EDE2ED1DEDE212E2121D12E2EDE212E2 

2  256  EDE2ED1DEDE212E2EDE2ED1D121DED1DEDE2ED1DEDE212E2121D12E2EDE212E2 

EDE2ED1DEDE212E2EDE2ED1D121DED1D121D12E2121DED1DEDE2ED1D121DED1D 

4  256  FAC9F5C6FA36F539FAC90A39FA360AC6FAC9F5C605C90AC6FAC90A3905C9F539 

FAC9F5C6FA36F5390536F5C605C9F539FAC9F5C605C90AC60536F5C6FA360AC6 
FAC9F5C6FA36F539FAC90A39FA360AC605360A39FA36F5390536F5C6FA360AC6 
FAC9F5C6FA36F5390536F5C605C9F53905360A39FA36F539FAC90A3905C9F539 
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III.D  Proofs  of  Theorems  2  and  3 


It  is  convenient  to  develop  some  intermediate  results  for  a  family  of  related  systems.  Let 
denote  a  linear  system  which  is  the  cascade  of  a  rate-M  upsampler,  an  LTI  system 
whose  unit-sample  response  is  /^^[^^],  a  linear  time  varying  system  whose  kernel  is  fij[n;k], 
another  LTI  system  with  unit-sample  response  hj[—n],  and  finally  a  rate-M  downsampler. 
Hence 


q[n]  —  Sij{p[n]}  —  ^  1]  hj\m  —  nM\  hi[m  —  I  —  kM\p{k],  (111.20) 

m,l,k 

and  it  is  straightforward  to  show  that  the  kernel  of  this  system  is 

Uij[n-,k]  =  E  fij  [nM  -I-  m;  kM  -f  Z]  hj [m]  hi[m  —  Z].  (111.21) 

m,l 

We  begin  with  the  following  lemma. 

Lemma  4  Let  hi[n]  and  hj[n]  be  chosen  from  an  set  of  M  orthogonal  signatures,  and 
let  /ij[n;A:]  for  j  =  1,2,  ...,M  be  an  admissibly  ergodic  family  of  kernels  in  the  sense  of 
Definition  2.  Then,  with  Vh  and  Xh  defined  in  (11.21)  and  (11.29),  respectively,  as 
Vfi  oo  and  Xh  oo,  the  kernel  (III. 21)  obeys 

Uij  [n;  k]^E  [Fii]  ^ [A:]  S[i  -  j]  =  Uy  [A:]  (111.22) 

and 

Uij[n;  n  —  A:]  u*jr  [m;  m  —  k]  ~  ~  (III.23) 

k  ^ 

where 

Uij[n-,  k]  =  Uij[n-,  k\  -  Uij[k]. 

Proof: 

Using  (II.44a)  with  (III.21),  we  obtain 

E  [uijim  A:]]  =  E  [F<j]  J]  hj[l]hi[l  -I-  kM] 

i 

and 

Uij[n-,  A:]  =  fij[nM  -I-  m;  kM  4-  Z]  hj[m]  hi]jn  -  Z] 

rriyl 

where  /ij[n;A:]  is  as  defined  in  (11.41). 

Prom  (III.25)  we  get 

E  [\uij[n-,  A:]|^]  =  ^2  ^ij  ~  i^'ykM  -  Z]  hj{m]  hj[m']  hi[l  -I-  m]  hi[l  -t-  m'], 

m,m\/ 


(in.24) 

(111.25) 

(111.26) 
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where  R^[n;  fe]  is  as  defined  in  (II.44b).  Applying,  in  order,  the  triangle  inequality  and  the  Cauchy 
inequality,  (111.26)  can  be  bounded  by 

E  WuiAn-,  jfclPl  <  ipihi)  (p(hj)  (III.27) 


where 


£  [|tty[n;jfc]p]  <  ip{hi)  (p{hj) 

(p‘^{g)=  \R^j[m-m']kM -l]\g'^[in]g^[rn']. 


(IIL28) 


Applying  the  Cauchy  inequality  again  to  ip^ig),  however,  gives,  after  some  simplification,  and  using 
(II.44d), 

v\9)<S^,dT^s 


where 


Hence. 


'^a= 


which  tends  to  zero  as  T>h  oo-  Collectively,  (III.24)  and  (III.29)  establish  (111.22). 
To  show  (III.23),  we  begin  by  noting  that 


(III.29) 


pii  [n,  m]  =  Y  [m-,rn-k] 

k 

=  Yj  Yh  +  s;  t]  fij,  [mM  +  s' ;  t']  hj  [s]  h'-  [s']  Qhi  [nM  +  s  —  t,  mM  +  s'  —  t'] 


SyS'  tyV 


where  0hi[n,m]  is  as  defined  in  (11.25). 
Using  (II.44b),  we  obtain 


E  [n,  m]|  =  Y^  [(^  ~  m)M  +  s  -  s';  t]  hj[s]  hj-  [s']  Qhi  ]nM  +  s-t,  mM  +  s'  -  t] 


which,  with 


p'ij  [n,  »7i]  =  +  s  -  s';  t]  hj[s]  hj'  [s']  Qhi  [nM  +  s-t,  mM  +  s'  -  t]  (III.30) 


and  using  (11.27),  can  be  rewritten  as 

E  [pif  [n,  m]]  =  pf  [n,  m]  +  j^Y  ^  H 

t  s 

=  p*j’  [n,  m]  +  j^var  [Fy]  S[j  -  j']  S[n  -  m].  (III.31) 

Now  applying  the  triangle  inequality  and  the  Cauchy  inequality  to  (III.31),  we  obtain  the  bound 

|py  [’^jHI  ^  (III.32) 
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where 


(ni.33) 


=  E  -m)M  +  s-  s';  t]  h,-[s]  hj. [s']| 

s,s' ,t 

and 

^  |/?*j  [(n  -  m)M  +  s  —  s';  f]  /ij[s]  hji  [s']  0|.  [nM  +  s  —  f,  mM  +  s'  —  t]. 

8,8',t 

Applying  the  Cauchy  inequality  to  (III.34)  yields 

?2  <  66 

where 

^3  =  ^  [(n  —  Tn)M  +  s  —  s';  t]  hj[s]  hji [s']  0^.  [nM  +  s  —  t,  mM  +  s'  —  t], 

s,s' ,t 

Similarly,  applying  the  Cauchy  inequality  to  (III.33)  yields 


<  [  E  +  s  -  s';  t]  h)[s]  j 

\s,s',t  ) 


Similarly,  applying  the  Cauchy  inequality  to  (III.36)  yields 


(111.34) 

(111.35) 

(ni.36) 

(III.37) 


^3  <  (  ^  ji?]j  [(n  -  m)M  +  s  —  s';  t]  h|[s]  0^.  [nM  +  s  —  t,  mM  +  s'  —  t]  j  . 

\s,s',t  J 

■  I  ^  j  [{n  -  m)M  +  s  -  s';  t]  h/j,  [s']  0^.  [nM  +  s  —  t,  mM  +  s'  —  t]  |  . 

\s,s',t  ) 


(ni.38) 


where  the  last  inequality  in  (III.38)  results  from  using,  in  order,  the  simple  bound 


tj 


and  (11.30). 

Using  (in.37),  (III.38)  and  (III.35)  with  (III.39)  we  get 


\pij 


(III.39) 


which  tends  to  zero  as  x/i  oo- 
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Next,  we  define 


Py  [n,  m]  =  p -j'  [n,  m]-E  [pV.  [n,  m] 


and  note 


(III.40) 


[n, m]  =  ^  ^  [s,  s';  k,  k']  Qh,  [«  -  k;  s'  -  A:']  hj[s  -  nM]  hf  [s'  -  mM] 

s,s'  k,k' 

where  [s,  s';t,t']  is  as  defined  in  (11.43)  and  (11.42). 

In  turn,  using  (II.44c)  with  (III.40),  we  get 


E  |pf[n,m]|  =  E.vE«.|;  Tif[s-s',t-f,s-t^k,k',l,l'] 

/ij[s  —  nM]  hf[s'  -  mM]  hj[t  —  nM]  hj'[t'  -  mM] .(III.41) 


Now,  from  a  simple  application  of  the  Cauchy  inequality  to  (11.25)  we  obtain 

10fcjTi,m]|  <  1. 


(III.42) 


Applying,  in  turn,  the  triangle  inequality,  the  bound  (III.42),  and  the  Cauchy  inequality  to  (III.41) 
we  obtain 

,2  (  ,  \ 

E  |p’{  [n,  m]  1^  <  XI ®  4  [«'  - 


X  X  [s-s',t-t',s-t-,  k,  k',  I,  I']  h) [t  -  nM]  h^,  [f'  -  mM] 


s,s' 

k,k'  1,1' 


Applying  the  Cauchy  inequality  once  again  to  the  right-hand  side  terms  of  (111.43)  yields 

■  . ./  21  ^ 

E  pljln,m]  < 

fxX  T^'ls-s',t-t’,s-t;k,k',l,l']l  h,^[s-nM]^ 


(III.43) 


(III.44) 


s.s'  i,t> 
k,k'  l,t' 


a, s' 

k,k'  1,1' 


( 

X  X  [s-  ,s-t\  k, k',l, Z']|  h|[t  -  nM] 

,  s.s'  t.t' 

\k.k'  l.l' 
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(III.45) 


Proposition  1  Suppose  that  the  input  to  the  system  defined  by  (III. 21)  is  a  zero-mean, 
white  Gaussian  sequence  Pi[n]  with  variance  cf,  and  let  the  corresponding  output  sequence 


be  qij[n].  Moreover,  assume  that  for  different  values  of  i  the  sequences  Pi[ri\  are  mutually 
independent.  Then  given  the  same  hypotheses  of  Lemma  4,  we  have 

qij  [n]  ^  E  [Pji]  <5 [i  -  j]  pi  [n]  +  Zij  [n]  (III.46) 

where  the  Zij\ri\  are  mutually  uncorrelated,  zero-mean  white  marginally  Gaussian  sequences 
with  variances  crfvar [F^j/M,  i.e., 

E[zij[n]]  =  0  (IIL47a) 

E[zij[n]pl.[m]]  =  0  (III.47b) 

E  [zij [n]  Ziijt  [m]]  =  -  i']  6\j  -  j']  (5[n  -  m]  -^aivai  [Fjy].  (III.47c) 

Furthermore,  we  have  that  the  sequences  9y[n]  are  zero-mean  and  white  with  variance 

var  qij [n]  =  -  jjcrf  \E  [Fj*]  p  +  ^crf  var  [Fy ]  (III.48) 


and  are  mutually  uncorrelated  for  different  values  of  j. 

Proof: 

We  begin  by  noting  that 

k 

can  be  rewritten  as 

qij[n]  =  E [Fii]  5[i  -  j]p[n]  +  Zifin] 

where 

^ij  W  =  Wy  [«;  h]  Pi  [n  -  fc] .  (III.49) 

k 

Prom  (III.49)  we  get  immediately,  since  zero-mean,  that  Zij[n]  satisfies  (III.47a),  i.e., 

Ep  ~  0) 

where  we  use  Ep  [•]  to  denote  expectation  with  respect  to  p[n]  given  fixed  but  arbitrary  realizations 
of  the  kernels  /ij[n;  k]. 
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In  addition,  from  (111.49)  we  get 

Ep  [zij  [n]  p*,  [m]]  =  cr?  S[i  -  i']  ti[n;  n  -  m].  (III.50) 

Hence,  since  (III.22)  in  Lemma  4  implies 

n[n;  k]  ^  0, 

we  get,  from  (III.50),  the  result  (in.47b),  i.e. 

Ep  [zij[n]pi,[m]]  ^  0. 

Finally,  using  (111.49)  we  also  obtain 

Ep  [2y [n]  z’j, [m]]  =  a?  5[i  -  i']  k] [m; l]Ep  |pi[n  -  k]pi[m  -  i]] 

k,l 

=  -  i']  ^ Uij[n-,  n-l]  ulj,  [m; m-l].  (111.51) 

k 

Applying  (III.23)  in  Lemma  4  to  (111.51)  we  obtain  (III.47c)  immediately,  i.e., 

Ep  [zij[n]  4,[m]]  c7?  S[i  -  *'];^var  [Fy]  5[n  -m]5\j-  j']. 

Finally,  to  show  that  Zij[n]  is  marginally  Gaussian  requires  a  straightforward  Central  Limit 
Theorem  argument. 


III.D.l  The  Forward  Link  Theorem 

Using  Proposition  1,  we  can  readily  establish  Theorem  2.  In  particular,  due  to  linearity 
we  can  partition  Xm[n]  into  two  components:  Xm^[n],  which  is  generated  by  the  set  of 
transmitted  sequences  Xj[n],  and  x^^[n],  which  is  generated  by  the  background  noise  it;[n]. 

We  first  note  that  if  in  Proposition  1  we  let  and  [n;  k]  =  c[n;  fc],  then  we 

readily  obtain,  using  superposition,  that 

M 

^  F [q Xm[n]  +  r;^)[n],  (III.52) 

i=l 


where 


M 


i=l 


Note,  in  addition,  that  due  to  the  properties  of  the  Zy  [n]  in  Proposition  1  the  are  mu¬ 

tually  uncorrelated,  zero-mean,  white  marginally-Gaussian  noise  sequences  with  variances 


,  M 

var  [n]  =  var  [C] 


2=1 


(III.53) 
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Next  we  note  that  if  in  Proposition  1  we  let  pi[n]  =  ^  u;[A:]  hm[k  -  nM]  and  [n;  k]  = 

k 

6[n;  fc],  then  we  again  readily  obtain  that 

M 

N  =  E  ^irn[n]  ^  [n],  (111.54) 

Z=1 

(2^ 

where  the  Vm  [n]  are  mutually  uncorrelated,  zero- mean,  white  marginally  Gaussian  noise 
sequences  with  variances 

W  =  A/bWoljB  [B]  p  -f- ./Vo Wovar  [B]  =  MqWqE  [|Bp]  .  (III.55) 

Hence,  combining  (III.52)  and  (III.54)  we  obtain  (11.46)  where 

^m[n]  =  v^^[n]-¥v^^[n]. 

Furthermore,  since  the  Xm[n]  and  w\p\  are  uncorrelated,  we  obtain  (11.47)  using 

var  UmN  =  var  [n]  +  var  [n] 

with  (111.53)  and  (III.55). 

III.D.2  The  Reverse  Link  Theorem 

Using  Proposition  1,  we  can  also  readily  establish  Theorem  3.  Again,  exploiting  hnearity  we 
partition  Xrn[n]  into  two  components:  x^\n],  which  is  generated  by  the  set  of  transmitted 
sequences  Xi[n],  and  x^[n],  which  is  generated  by  the  background  noise  zu[n]. 

When  in  Proposition  1  we  let  Pi[n]  =  Xi[n]  and  fij[n-,k]  =  Cij[n\k],  we  readily  obtain, 
using  superposition,  that 

M 

W  =  E  W  ^  ^  [Cmm]  Xm  N  +  [«] ,  (III.56) 

2=1 

where 

M 

=  E  ^irn[n]. 

i=l 

Note,  in  addition,  that  due  to  the  properties  of  the  Zij[n]  in  Proposition  1,  the  are 

mutually  uncorrelated,  zero-mean,  white  marginally- Gaussian  noise  sequences  with  vari¬ 
ances 

j  M 

vax  [n]  =  —  E  var  [Cim].  (III.57) 

2=1 

Next,  when  in  Proposition  1  we  let  pi[n]  =  E^W  ^rn{k  -  nM]  and  fij[n-,k]  =  bj[n;k], 

k 
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we  again  readily  obtain  that 


M 

W  =  M  ^  [n] ,  (111.58) 

i—\ 

where  the  v^[n\  are  mutually  uncorrelated,  zero-mean,  white  marginally  Gaussian  noise 
sequences  with  variances 

var [n]  =  MoWq\E  [B,^]  +  A^oWovar  [Bm]  =  NqWqE  [\Bm\^] .  (111.59) 

Hence,  combining  (III.56)  and  (III.58)  we  obtain  (11.55  where 

Vm[n]  =  +  u^)[n]. 

Furthermore,  since  the  Xm{^  and  w[n]  are  uncorrelated,  we  obtain  (11.56)  using 

var  Vm[n]  =  var  vg}  [n]  -I-  var vg}  [n] 

with  (III.57)  and  (III.59). 


(111.60) 

(111.61) 

(111.62) 

(111.63) 

(111.64) 


Note  that  by  symmetry,  Ai,A2,... ,  Am  are  also  zero-mean,  mutually  uncorrelated  random 
variables,  whose  variances  we  denote  by  Ai,  A2, . . . ,  \m- 
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Now,  any  B  with  finite  variance  can  be  expanded  in  the  form 


B  =  £  +  Y,  VkAk 


(111.65) 


where,  for  all  k,  E  [£*Afc]  =  0  and  %  are  complex  constants.  In  particular,  it  suffices  to 
choose 

r)k  =  E[B*Ak]  IE  [lifcpj. 

Furthermore,  from  (11.58)  we  see  that  if  B  maximizes  7,71,  so  does  kB  for  any  k.  Hence,  to 
fix  a  particular  solution,  we  may,  without  loss  of  generality,  set  77^  =  1. 

Using  (III.65),  (III.64)  simplifies  to 


7m(5)  =  XI,{e  [\£f]  +  Xm{l  -  Xm/M)  +  Y,  (1  -  Xk/M) . 


which,  since 


Xi  =  E 


is  maximized  when  e  =  0  and  %  =  0  for  A:  7>^  m.  Thus,  the  maximum  value  of  7^  is  obtained 
when  B  is  of  the  form  (11.59). 

To  obtain  the  bound  (11.61),  it  suffices  to  substitute  the  optimum  value  of  b  into  (III.60), 
which  yields 


where 


^  IM  =  _ 1 

^m  =  E[Sra\Am?ie] 


It  therefore  remains  only  to  note  that 


■A/b>Vo  +  j^'Y  ^k\-Ak\^ 
l-cpm/M  =  E  - - 

^foWo  +  JJY^k\Ak\^ 


III.F  Proof  of  Lemma  2 

Let  ^fn  denote  the  left  hand  side  of  (11.67).  Then 


=  6  =  •  •  •  =  Cm  =  ^0 
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and,  hence, 


««  - 

=  —E 
M 


L  k 

M  +  (M-l)Efc^ifc 


1  +  Efc  t’fc 


M-l  1 


+  77^’ 


where 


with 


M  M 


(p  =  E 


1+7, 


7  = 

k 

Now  7  as  defined  in  (III.68)  is  an  Erlang  random  variable  of  order  M  and 
Hence, 


-I 


00  ^M;yM-lg-f£7 

(1  +  7)  (M-l)! 


dj. 


Exploiting  the  identities  [15] 


fOO  p  —  SZ 

rOO 

Jo 


and 


we  get 


fc— 1 

+ 

fc=0 


4.  V' /'_l  'sn-k-l  fk 

t  +  1  t+1  ‘ 


p  = 


M 


M 


M-2 


fc:=0 


(M  - 1)! 

Finally,  substituting  (III. 70)  into  (III.66)  we  get  (11.67). 


(111.66) 

(111.67) 

(111.68) 
mean  M/)U. 

(111.69) 


(III.70) 
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